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ABSTRACT 

We present a new method to achieve high-contrast images using segmented and/or on- axis tele- 
scopes. Our approach relies on using two sequential Deformable Mirrors to compensate for the large 
amplitude excursions in the telescope aperture due to secondary support structures and/or segment 
gaps. In this configuration the parameter landscape of Deformable Mirror Surfaces that yield high 
contrast Point Spread Functions is not linear, and non-linear methods are needed to find the true 
minimum in the optimization topology. We solve the highly non-linear Monge- Ampere equation that 
is the fundamental equation describing the physics of phase induced amplitude modulation. We de- 
termine the optimum configuration for our two sequential Deformable Mirror system and show that 
high-throughput and high contrast solutions can be achieved using realistic surface deformations that 
are accessible using existing technologies. We name this process Active Compensation of Aperture 
Discontinuities (ACAD). We show that for geometries similar to JWST, ACAD can attain at least 
10~ ' in contrast and an order of magnitude higher for both the future Extremely Large Telescopes and 
on-axis architectures reminiscent of HST. We show that the converging non-linear mappings resulting 
from our Deformable Mirror shapes actually damp near-field diffraction artifacts in the vicinity of the 
discontinuities. Thus ACAD actually lowers the chromatic ringing due to diffraction by segment gaps 
and strut's while not amplifying the diffraction at the aperture edges beyond the Fresnel regime. This 
outer Fresnel ringing can be mitigated by properly designing the optical system. Consequently, ACAD 
is a true broadband solution to the problem of high-contrast imaging with segmented and/or on-axis 
apertures. We finally show that once the non-linear solution is found, fine tuning with linear methods 
used in wavefront control can be applied to further contrast by another order of magnitude. Generally 
speaking, the ACAD technique can be used to significantly improve a broad class of telescope designs 
for a variety of problems. 

Subject headings: planetary systems - techniques: coronagraphy, wavefront control 



1. INTRODUCTION 



Exo-planetary s ystems that are directly imaged usin. 



existing facilities (Marois et al. 2008 Kalas et al.||2008 



ng 
38; 



Lagrange et al. 201 0| ) give a unique laboratory to con 



strain planetary fo rmation at wide separations (Rafikov 
20051 IDodson-Robinson eTal^|20l)9l |Kratter et al.p010| 



Johnson et al.||2010p , to study the p lanetary luminosity 
distribution at critical yo ung ages ( jSpiegel fc Burrows 



2012 Fortney et al. 2008) and the atmospheric proper 



ties ofTow surface "gravity objects ( Barma n et al. 2011b|a 



Marley et aTpOloj [Madhusudhan et al.||201ip . Upcom 



ing surveys, conducted wi th instruments specifically de- 
signed for high-contrast (Dohlen et al. 2006[ | Graham] 



et al. 2007 Hinkley et al. '2TjTT), will unravel the bulk 



of this population of self-lummous jovian planets and 
provide an unprecedented understanding of their forma- 
tion history. Such instruments will reach the contrast 
required to achieve their scientific goals by combining 
Extreme Ad aptive Optics systems (Ex- AO, |Poyneer fc 



Veran (2005)), optimized coronagraphs ( |Soummer et al 



2011| |Guyon||2003| |Rouan et a l.||2000 [) and nanometer 



class w avefront calibration (Sauvage et al. 2007 Wal 



lace et al.|[2009| |Pueyo et al.|]2010[ ). In the future higE 



contrast instruments on Extremely Large Telescopes will 
focus on probing planetary formation in distant star 
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forming regions (Macintosh et al. 2006 ), characterizing 
both the spectra ot cooler gas giants (Verinaud et al 



2010 1 and the reflected light of planets in the habitable 
zone of low mass stars. The formidable contrast neces- 
sary to investigate the presence of biomarkers at the sur- 
face of earth analogs (> 10 10 ) cannot be achieved from 
the ground beneath atmospheric turbule nce and will r e- 
quire dedicated space-based instruments (Guyon||2005). 



The coronagraphs that will equip upcoming Ex- AO in- 
struments on 8 meter class telescopes have been designed 
for contrasts of at most ~ 10~ 7 . Secondary support 
structures (or spiders: 4 struts each 1 cm wide, ~ 0.3% 
of the total pupil diameter in the case of Gemini South) 
have a small impact on starlight extinction at such lev- 
els of contrasts. In this case, coronagraphs have thus 
been optimized on circularly symmetric apertu res, which 
only take in to account the central obscuration ( Soummer 
et al.poTT l. However, high-contrast instrumentation on 



tuture observatories will not benefit from such gentle cir- 
cumstances. ELTs will have to support a substantially 
heavier secondary than 8 meter class observatories do, 
and over larger lengths: as a consequence the relative 
area covered by the secondary support will increase by a 
factor of 10 (30 cm wide spiders, occupying ~ 3% of the 
pupil diameter in the case of TMT) . This will degrade the 
contrast of coronagraphs only designed for circularly ob- 
scured geometries by a factor ~ 100, when the actual en- 
visioned contrast for an ELT exo-planet imager can be as 
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low as ~ 10" 8 ( [Macintosh et al.|2006| ). While the trade- 
offs associated with minimization ot spider width in the 
space-based case have yet to be explored, secondary sup- 
port structures will certainly hamper the contrast depth 
of coronagraphic instruments of such observatories at lev 



jugate to the telescope pupil (Shaklan & Green 2006 
Pueyo & Kasdin| |2007|) . Indeed a single DM can only 



mimic halt ot the spatial frequency content of amplitude 
errors and compensate for them only on one half of the 
image plane (thus limiting the scientific field of view) 
over a moderate bandwidth. In theory, architectures 
with two sequential Dcformable Mirrors, can circumvent 
this problem and create a symmetric broadband high 
contr ast PSF ( |Shaklan fc Green||2006[ |Pueyo fc Kasdin 
( Guyon et al. 2008 Trauger et al.|[20i0 1 . Coronagraphs 2007 1 . The first demonstration of s ymmetric dark hole 



els that are well above the 10 contrast requirement. 
As a consequence, telescope architectures currently en- 
visioned for direct characterization of exo-earths consist 
of monolithic, off-axis, and thus un-obsc ured, telescopes 



for such architectures take advantage ot the pupil symme- was reported in Pueyo et al.| J2009j) and has s ince been 



try t o reach a theoretical contrast of ten orders of magni- 



tude (|Guyonetal.|2005 


Vanderbei et al. 2003a b Kasdin 


et al. ||2005t IMawet et a 


. kuchner & Tr'auk 2002; 


Soummer et al.||2003|. J 


however, using obscured on- axis 



limited real estate associated with a given launch vehicle 
and can allow larger apertures that increase the scien- 
tific return of space-based direct imaging survey. Recent 
solutions can mitigate the presence of secondary support 
structures in on-axis apertures. However these concepts 
prese nt practical limitation s: APLCs on arbitr ary aper 



tures ( Soummer et al.|2009| and Shaped Pupils ( Carlotti 
201 1| ) sutler from throughput loss for ve ry hign 



et al. 

contrasTcIesigns, and PIAAMCM (|Guyon et aL||2010al 



rely on a phase mask technology whose chromatic prop 
erties have not yet been fully characterized. Moreover 
segmentation will further complicate the structure of the 
telescope's pupil: both the amplitude discontinuities cre- 
ated by the segments gaps and the phase discontinuities 
resulting from imperfect phasing will thus further de- 
grade coronagraphic contrast. Devising a practical so- 
lution for broadband coronagraphy on asymmetric, un- 
friendly apertures is an outstanding problem in high con- 
trast instrumentation. The purpose of the present paper 
is to introduce a family of practical solutions to this prob- 
lem. As their ultimate performances depend strongly on 
the pupil structure we limit the scope of this paper to 
a few characteristic examples. Full optimization for spe- 
cific telescope geometries can be conducted as needed. 

The method proposed here takes advantage of state- 
of-the art Dcformable Mirrors in modern high-contrast 
instruments to address the problem of pupil amplitude 
discontinuities for on-axis and/or segmented telescopes. 
Indeed, coronagraphs are not sufficient to reach the high 
contrast required to image faint exo-planets: wavefront 
control is needed to remove the light scattered by small 



impe rfections on the optical surfaces (Brown & Burrows 
19901. Over the past eight years, significant progress has 



been made i n this area, both in the development of new 
algorithms ( |Borde fc Traub||2006| |Give'on et"aL[|2007 l 
and in the experimental demonstration o t high-contrast 
imaging with a variety of coronagraphs flGive'on et ah] 
"20071 [Trauger fc Traub|2007]|Guyon et al.|20l0b[|Belikoy| 
et al.||201ip . These experiments rely on a system with 
a single Deformable Mirror which is controlled based on 
diagnostics downstream of the coronagraph, either at the 
science camera or as close as possible to the end detector 



(Wallace et al. 
rations are wei 



2009 Pueyo et al. 20101. Such configu- 



suited to correct phase wavefront errors 



arising from surface roughness but have limitations in 
the presence of pure amplitude errors (reflectivity), or 
phase-induced amplitude errors, which result from the 
propagation of surface errors in optics that are not con- 



generalized to broadband by |Grofl et al. (2011). In such 
experiments the coronagraph has been designed over a 
full circular aperture, the DM control strategy is based 
on a linearization of the relationship between surface de- 
formations and electrical field at the science camera, and 
the modeling tools underlying the control loop consist 
of classic Fourier and fresnel propagators. This is illus- 
trated on the left panel of Fig. [I] As a consequence, a 
wavefront control system composed of two sequential De- 
formable Mirrors is currently the baseline architecture of 
current ly envisioned coronagraphic space-bas ed instru- 
ments ( [Shaklan et "aT1|2 006 Kri st et al |[20lT] ) and ELT 
planet imagers ( |Macmtosh et aiT||200fc>J ). Une can thus 
naturally be motivated to investigate if such wavefront 
control systems can be used to cancel the light diffracted 
by secondary supports and segments in large telescopes, 
since such structures are amplitude errors, albeit large 
amplitude errors. 

The purpose of our study is to demonstrate that in- 
deed a two Deformable Mirror (DM for the remainder 
of this paper) wavefront control system can mitigate the 
impact of the pupil asymmetries, such as spiders and 
segments, on contrast and thus enable high contrast on 
unfriendly apertures. In §. [2] we first present a new ap- 
proach to coronagraph design in the presence of a central 
obscuration, but in the absence of spiders or segments. 
We show that for coronagraphs with a pupil apodization 
and an opaque focal plane stop, contrasts of 10 -10 can be 
reached for any central obscuration diameter, provided 
that the Inner Working Angle is large enough. Natu- 
rally the secondary support structures, and in the seg- 
mented cases, segment gaps, will degrade this contrast. 
As our goal is to use two DMs as an amplitude modu- 
lation device, we first briefly review in §. [3] the physics 
of such a modulation. In §. [4] we introduce a solution 
to this problem: we show how to compute DM surfaces 
that mitigate spiders and segment gaps. Current algo- 
rithms used for amplitude control operate under the as- 
sumption that amplitude errors are small, and thus they 
cannot be readily applied to the problem of compensat- 
ing aperture discontinuities, which have inherently large 
reflectivity non-uniformities. Fig. [I] illustrates how the 
present manuscript introduces a control strategy for the 
DMs that is radically different from previously published 
amplitude modulators in the high-contrast imaging lit- 
erature. Our technique, which we name Active Com- 
pensation of Aperture Discontinuities (hereafter ACAD), 
finds the adequate DM shapes in the true non-linear large 
amplitude error regime. In this case the DMs' surfaces 
are calculated as the solution of a non-linear partial dif- 
ferential equation, called the Monge-Ampere Equation. 
We describe our methodology to solve this equation in 
§. [4] and illustrate each step using an obscured and seg- 
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mented geometry similar to JWST. As ACAD DM sur- 
faces are prescribed in the ray-optics approximation this 
is a fundamentally broadband technique provided that 
chromatic diffractive artifacts, edge ringing in particu- 
lar, do not significantly impact the contrast. This is 
what we discuss in § . [5] We find that when remapping 
small discontinuities with Dcformablc Mirrors, the spec- 
tral bandwidth is only limited by wavelength-dependent 
edge-diffrac tion ringing in the Fresn el approximation (as 
discussed in Pueyo & Kasdin ( 2007[ ) for instance). High- 
contrast instruments where this ringing is mitigated have 
already been designed; while future work in high preci- 
sion optical modeling is necessary to fully quantify the 
true chromatic performances of ACAD, we do not expect 
these effects to be a major limitation to broadband oper- 
ations. In §.[6] we present the application of our method 
to various observatory architectures. Note that the con- 
trast levels stated in §. [6] represent a non-optimal esti- 
mate of ACAD performances with on-axis and/or seg- 
mented apertures. Our calculations are carried out in 
the absence of atmospheric turbulence, quasi-static wave- 
front errors or coronagraphic manufacturing defects. We 
discuss these limitations in §. [7J with a specific empha- 
sis on quasi-static phase errors in a segmented telescope. 
We show that when the aperture discontinuities are thin 
enough field distortion is negligible for spatial frequen- 
cies within the field of view defined by the DMs control- 
lable spatial frequencies. We then discuss issues associ- 
ated with phase discontinuities when applying ACAD to 
a segmented telescope. We show that they can be cor- 
rected by superposing single DM classical wavefront con- 
trol solutions to the ACAD shape of the second DM. We 
finally argue, that should high precision diffractive mod- 
els be developed, then the solutions presented herein can 
be used as the starting point of dual DMs iterative algo- 
rithms relying on an image-plane based metric and thus 
lead to higher contrast than reported herein. Most of the 
future exo-planet imagers, either on ELTs or on future 
space missions, are envisioned to control their wavefront 
in real time with two sequential DMs. The method pre- 
sented in this manuscript thus renders high contrast coro- 
nagraphy possible on any observatory geometry without 
adding any new hardware. 

2. CORONAGRAPHY WITH A CENTRAL OBSCURATION 

2.1. Optimizing pupil apodization in the presence of a 
central obscuration 

Because the pupil obscuration in an on-axis telescope 
is large it will be very difficult to mitigate its impact with 
DMs with a limited stroke. Indeed, the main hindrance 
to high-contrast coronagraphy in on-axis telescope is the 
presence of the central obscuration: it often shadows 
much more than 10% of the aperture width while sec- 
ondary supports and segments gaps cover ~ 1%. We 
thus first focus of azimuthally symmetric coronagraphic 
designs in the presence of a central obscuration. This 
problem (without the support structures) has been ad- 
dressed in previous publicatio ns either using circula rly 
symmetric pupil apodiz ation (jSoummer et al. 2011 ) or 
a series of phase masks (Mawet et al.||20Il|). Both solu- 



tions however are subject to limitations. The singularity 
at the center of the Optical Vector Vortex Coronagraph 
might be difficult to manufacture and a circular opaque 
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Figure 1. Top, blue: Envisioned architecture of future exo-earth 
imaging missions: a monolithic un-obscured telescope feeds a coro- 
nagraph designed on a circular aperture. The wavefront errors are 
corrected using two sequential Deformable Mirrors (DMs) that are 
controlled using a quasi-linear feedback loop based on image plane 
diagnostics. The propagation between optical surfaces, between 
the DMs in particular is assumed to occur in the Fresnel regime. 
Bottom, orange: ACAD solution: the coronagraph is designed 
for a circular geometry around the central obscuration. The two 
sequential DMs are controlled in the non-linear regime based on 
a pupil plane cost function. The propagation between the DMs is 
modeled using ray optics, and a we conduct a quantitative one- 
dimensional analysis of the diffraction artifacts. 

spot thus lies in the central portion of the phase mask 
(< X/D) which results in a de gradation of the ide al con- 



trast of such a coronagraph (Krist et al. 2011). The 



solutions in Soummer et al. (2011) result from an opti 



mization seeking to maximize the off-axis throughput for 
a given focal plane stop diameter: the final contrast is 
absent from the optimization metric and is only a by- 
product of the chosen geometry. Higher contrasts are 
then obtained by increasing the size of the focal plane 
mask, and thus result in a loss in IWA. 

2.2. The optimization problem 
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Here we revisit the solution proposed by Soummer 
et al. (2011 1 in a slightly different framework. We recog- 



mze that, m the presence of wavefront errors, high con- 
trast can only be achieved in an area of the field of view 
that is bounded by the spatial frequency corresponding 
to the DM's actuator spacing. We thus consider the 
design of an Apodized Pupil Lyot Coronagraph which 
only aims at generating high contrast between the Inner 
Working Angle (IWA) and Outer Working Angle (OWA). 
In order to do so , we r ewrite coronagraphs described by 
Soummer et al. (2003) as an operator C which relates 
the entrance pupil ~PJr) to the electrical field in the final 
image plane. We first call the Hankel transform of 
the entrance pupil: 



D/2 



Ds/2 



P(r)J {r£)rdr 



where D is the pupil diameter, D$ the diameter of the 
secondary and £ the coordinate at the science detector 
expressed in units of angular resolution (Xq/D). Ao is the 
design wavelength of the coronagraph chosen to translate 
the actual physical size of the focal plane mask in units of 
angular resolution (often at the center of the bandwidth 
of interest). A is then the wavelength at which the coro- 
nagraph is operating (e.g. the physical size of the focal 
plane mask remains constant as the width of the diffrac- 
tion pattern changes with wavelength). For the purpose 
of the monochromatic designs presented herein A = Ao . 
Then the operator is given by: 



C[P(r)] (0 = P(0~ 



A^ 
A?, 



D/2 



P(ri)K(S,Ti)TidT) (2) 



where K(£, jf) is the convolution kernel that captures the 
effect of the focal plane stop of diameter M stop - 



»/ 2 



Jq(ut]) Jo(u^)udu 



(3) 



An analytical closed form for this kernel can be calcu- 
lated using Lommel functions. Note that this Eq. [2] as- 
sumes that the Lyot stop is not undersized. Since we are 
interested in high contrast regions that only span radially 
all the way up to a finite OWA, we seek pupil apodization 
of the form: 



P(r)= J2 PM^) 

i n U 



(4) 



fc=0 



where Jq{t) denotes the Bessel function of the first 
kind of order Q and a? the feth zero of this Bessel 
function. In order to devise optimal apodizations over- 
obscured pupils, Q, can be chosen to be large enough 
so that Jq(t) <C 1 for r < Ds/2 (in practice we choose 

Q = 10). The a J? corresponds to the spatial scale of 
oscillations in the coronagraph entrance pupil, and such 
a basis set yields high contrast regions all the way to 
OWA ~ N mo( ies^o/D. Since the operator in Eq. [2] is 
linear, finding the optimal pt can be written as the fol- 



lowing linear programming problem: 

max min(P(r)) Under the constraints: (5a) 
{Pk} L r 1 

\C [{P{r))\ (01 < 1(rv ^ for IWA < £ < OWA ( 5b ) 

max(P(r)) = 1 (5c) 

r 

, d 



dr 



[P(r)]\<b. 



(5d) 



Our choice of cost function and constraints has been 
directed by the following rationale: 

5. a We maximize the smallest value on the apodization 
function in an attempt to maximize throughput. 
The actual throughput is a quadratic function of 
the pk- Maximizing it requires the solution of a 
non-linear optimization problem (as described in 



Vanderbei et al.| ( |2003b| >) 



5.b The contrast constraint is enforced between the 
IWA and the OWA (< N modes ). 

5.c The maximum of the apodization function is set 
to one (otherwise the pk will be chosen to be suffi- 
ciently small so that the contrast constraint is met). 

5.d The absolute value of the derivative across the 
pupil cannot be larger than a limit, denoted as b 
here. As the natural solutions of such problem are 
very osci llatory (or "bang bang" , |Vanderbei et al. 
(2003b]a)), a smoothness constraint has to be en- 
forced (see Vanderbei et al.] ([2007 ) for a similar 
case) . 

Note that the linear transfer function in Eq. [2] can also 
be derived for other coronagraphs, with grayscale and 
phase image- plane masks, or for the case of under-sized 
Lyot stops. As general coronagraphic design in obscured 
circular geometries is not our main purpose, we limit the 
scope of the paper to coronagraphs represented by Eq. [2] 

2.3. Results of the optimization 

Ty pica l results of the monochromatic optimization in 
Eqs. |5d| with A = Ao, arc shown in Fig. [2] for central 
obscurations of 10, 20 and 30%. In the first two cases 
the size of the focal plane stop is equal to 3X/D, the 
IWA is 4A/L> and the OWA is 30X/D. As the size of 
the central obscuration increases the resulting optimal 
apodization becomes more oscillatory and the contrast 
constraint has to be loosened in order for the linear pro- 
gramming optimizer to converge to a smoother solution. 
Alternatively increasing the size of the focal plane stop 
yield smooth apodizers with high contrast, at cost in an- 
gular resolutions (bottom panel with a central obscura- 
tion of 30%, a focal plane mask of radius 4A/£>, an IWA 
of bX/D and an OWA of 30A/D). T hese trade-offs were 
described in Soummer et al.| ( 2011[ ), however our linear 
programming approach to the design of pupil apodiza- 
tions now imposes the final contrast instead of having 
it be a by product of fixed central obscuration and fo- 
cal plane stop. These apodizations can either be gener- 
ated using grayscale screen (at a cost in throughput and 
angular resolution) or a series of two aspherical PIAA 
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A/D in case of a PIAAC 




r/D 




A/D in case of an APLC 



Figure 2. Optimal design Apodized Pupil Lyot Coronagraphs on circularly obscured apertures. With fixed obscuration ratio, size of 
opaque focal plane mask, IWA and OWA, our linear programming approach yields solutions with theoretical contrast below 10 10 . All PSFs 
shown on this figure are monochromatic for A = Ao- When all other quantities remain equal and the central obscuration ratio increases (from 
10% in the top panel to 20% in the middle panel), then the solution becomes more oscillatory (e.g less feasible) and the contrast constraint 
has to be relaxed. Eventually the optimizer does not find a solution and the size of the opaque focal plane mask (and thus the IWA) has 
to be increased (central obscuration of 30% in the bottom panel). On the right hand side, we present our results in two configurations: 
when the apodization is achieved using at grayscale screen (APLC), "on-sky" A/D bottom x-axis, and when the apodization is achieved 
via two pupil remapping mirrors (PIAAC), "on-sky" X/D top x-axis. We adopt this presentation to show that ACAD is "coronagraph 
independent" and that it can be applied to coronagraphs with high throughput and small IWA. 
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mirrors (for better throughput and angular resolution). 
In order not to lose generality, we present our results on 
Fig. [^considering the two types of practical implementa- 
tions (classical apodization and PIAA apodization). In 
the case of a grayscale amplitude screen the angular res- 
olution units are as defined in Eq. [2] and the through- 
put is smaller than unity. In the case of PIAA apodisa- 
tion the throughput is unity and the angular resolution 
units have been magnified by the field inde pendent cen- 
troid based angular magnification defined in |Pueyo et al.| 
(2011 ). We adopt this presentation for the remainder of 
the paper where one dimensional PSFs will be presented 
with "APLC angular resolution units" in the bottom hor- 
izontal axis and'TIAAC angular resolution units" in the 
top horizontal axis. 

Note that this linear programming approach only op- 
timizes the contrast for a given wavelength. However, 
since the solutions presented in Fig. [2] feature contrasts 
below 10 10 , we choose not to focus on coronagraph chro- 
matic optimizations. Instead, in order to account for the 
chromatic behavior of the coronagraph, the monochro- 
matic simulations in § . [6] are carried out under the con- 
servative assumption that the physical size of the focal 
plane stop is somewhat smaller than optimal (or that the 
operating wavelength of the coronagraph is slightly off, 
A = 1.2Ao). As a consequence the raw contrast of the 
coronagraphs presented in §.j6jis ~ 10 -9 . Note that this 
choice is not representative ofall possible Apodized Pupil 
Coronagraph chromatic configurations. It is merely a 
shortcut we us e to cover the variety of cases presented in 
§. [6] In §. |7.2| we present a set of broadband simulations 
that include wavefront errors and the true coronagraphic 
chromaticity for a specific configuration and show that 
bandwidth is more likely to be limited by the spectral 
bandwidth of the wavefront control system than by the 
coronagraph. However, future studies aimed at defining 
the true contrast limits of a given telescope geometry will 
have to rely on solutions of the linear problem in Eqs. |5d| 
which has been augmented to accommodate for broach 
band observations. In theory, the method presented here 
can also be applied to asymmetric pupils. However, the 
optimization quickly becomes computationally intensive 
as the dimensionality of the linear programming increases 
(in particular when the smoothness constraint and the 
bounds on the apodization have to be enforced at all 
points of a two dimensional array) . This problem can be 
somewhat mitigate d when seeking for b inary apodiza- 
tions, as shown in Carlotti et al.| ( 2011 ), at a cost in 
throughput and angular resolution. 

3. PHYSICS OF AMPLITUDE MODULATION 

3.1. General equations 

We have shown in § . [2] that by considering the design 
of pupil apodized coronagraphs in the presence of a cir- 
cular central obscuration as a linear optimization prob- 
lem, high contrast can be reached provided that the focal 
plane mask is large enough. In practice, the secondary 
support structures and the other asymmetric discontinu- 
ities in the telescope aperture (such as segment gaps) will 
prevent such levels of starlight suppression. We demon- 
strate that well controlled DMs can circumvent the ob- 
stacle of spiders and segment gaps. In this section, we 
first set-up our notations and review the physics of phase 



to amplitude modulation. We consider the system rep- 
resented on Fig. [3] where two sequential DMs are located 
between the telescope aperture and the entrance pupil 
of the coronagraph. In this configuration, the telescope 
aperture and the pupil apodizer are not in conjugate 
planes. This will have an impact on the chromaticity of 
the system and is discussed in § . [5] Without loss of gener- 
ality we work under the "foldecr assumption illustrated 
on Fig. [5] where the DMs are not tilted with respect to 
the optical axis and can be consid ered as lenses of in- 
dex of refraction —1 (as discussed in |Vanderbei &: Traub 
(2005) . In the scalar approximation the relationship be- 
tween the incoming field, EjjMl(x,y), and the outgoing 
field, Edm2{x2,U2), is given by the diffractive Huygens 
Integral: 



E D M2(x2,1/2) 



iXZ 



E DM i{x,y)e i2 -? Q ( x <y< x ^dxdy 

(6) 

where A corresponds to the telescope aperture and 
Q(x, y, x 2l 2/2) stands for the optical path length between 
any two points at DM1 and DM2: 

Q(x,y,x 2 ,y 2 ) = h 1 (x,y)+S(x,y,x 2 ,y 2 )-h 2 (x 2 ,y 2 ) (7) 

S(x,y,x 2 ,y 2 ) is the free space propagation between the 
DMs: 



y/{x- x 2 ) 2 + {y - y 2 ) 2 + {Z + h x (x,y) - h 2 {x 2 ,y 2 )) 2 

(8) 

where Z is the distance between between the two DMs, 
h\ and h 2 are the shapes of DM1 and DM2 respectively 
(as shown on Fig. [5| and A is the wavelength. We rec- 
ognize that two sequential DMs act as a pupil remap- 



ping unit similar to PIAA coronagraph (|Guyon| 120031) 

i n n . ' T 1 J I I l' 



whose ray optics eq uations were first derived by |lraub" 



& Vanderbei ( 2003[ ). We briefly state the notatio n use' 
to describe s uch an optical system as introduced in |Pueyo 



et al. (2011): 



• For a given location at DM2, {x 2 ,y 2 ), the location 
at DM1 of the incident ray according to ray optics 
is given by: 



xi{x 2 ,y 2 ) = fi(x 2 ,y 2 ) 
yi(x 2 ,y 2 ) = gi{x 2 ,y 2 ). 



(9a) 
(9b) 



• Conversely, for a given location at DM1, (xi,yx), 
the location at DM2 of the outgoing ray according 
to geometric optics is given by: 



x 2 {xi,yi) = f2{xi,yi) 
2/2(^1,2/1) = g 2 (xx,yx). 



(10a) 
(10b) 



• Fermat's principle dictates the following relation- 
ships between the remapping functions and the 
shape of DM1: 



dhi 
dx 

dy 



_ xi - f 2 {xi,y 1 ) 
(xi,yi) Z 

= y\ - ggCgbj/i) 



(lla) 
(lib) 
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• Conversely, if we choose the surface of DM2 to en- 
sure that the outgoing on-axis wavefront is flat, we 
then find: 



dh2 
dx 
dh 2 
dy 



_ x 2 - fi(x 2 ,y 2 ) 
02,1/2) Z 

= V2 - 9l{x 2 ,y 2 ) 
02,1/2) z 



(12a) 
(12b) 



3.2. Fresnel approximation and Talbot Imaging 



In 



Pueyo et al. (2011) we showed that one could ap- 



proximate the propagation integral in Eq. [6] by taking in 
a second order Taylor expansion of Q(x, y, x 2 ,y 2 ) around 
the rays that trace (fi(x 2 , y 2 ),g , i(x 2 , y 2 )) to (x 2 ,y 2 ). In 
this case the relationship between the fields at DM1 and 
DM2 is: 



E 



DM2 



(x 2 ,y 2 ) 



e >■ 



iXZ 



E 



DM1 



(x,y)e' 



l( x -x 1 ) 2 +2'^(x-x 1 )(y-y 1 ) + i^(y^y 1 ) 



dxdy 



(13) 



When the mirror's deformations are very small compared to both the wavelength and D 2 /Z, the net effect of the 
wavefront disturbance created by DM1 can be captured in Edmi( x : v) an d the surface of DM2 can be factored out of 

dx \x-l,vi 



Eq. j6j In this case x\ — x 2 , y\ = y 2 

EdM2(X2,V2) 



1 dh 
1 ' dy 



= 0, 



= 1. Then, Eq. 



13 



e >■ 



{Z-h 2 (x 2 ,y 2 )) 



iXZ 



e >• 



CuO,2/)) fk((x-X2) 2 + (y-y 2 ) 2 ) 



reduces to: 
dxdy 



(14) 



Aperture 



which is the Fresnel approximation. If moreover 
hi(x,y) = Aecos(^j(rair + ny)), h 2 (x,y) = -hx(x,y), 
with e <C 1, then the outgoing field is to first order: 



E D M2(x2,y 2 ) oc 



n\Z(m 2 - 



n2 ) x / 27r / \ \ 
Aecos( — (mx 2 + ny 2 )). 

(15) 



This phase-to-amplitude coupling is a well known optical 
phenomenon called Talbot imaging and was introduced 
to the conte xt of high contrast imaging by Sha klan fc| 
Green ( 2006 ). In the small deformation regime, the phase 
on DM1 becomes an amplitude at DM2 according to the 
coupling in Eq. 15 When two sequential D Ms are con- 
trolled to cancel small amplitude errors, as in |Pueyo et al 



(2009), they operate in this regime. Note, however, that 



the coupling factor scales with wavelength (the resulting 
amplitude modulation is wavelength independent, but 
the coupling scales as A): this formalism is thus not ap- 
plicable to our case, for which we are seeking to correct 
large amplitude errors (secondary support structures and 



segments) with the DMs. In practice, when using Eq. 15 
in the wavefront control scheme outlined in |Pueyo et al 



(2009) to correct aperture discontinuities, this weak cou- 



pling results in large mirror shapes that lie beyond the 
range of the linear assumption made by the DM con- 
trol algorithm. For this reason, methods outlined on the 
left panel of Fig. [l] to correct for aperture discontinu- 
ities do not converge to high contrast. Because phase to 
amplitude conversion is fundamentally a very non -linear 
phenomena, these descending gradient methods ( Borde 



Bo 

fc Traub|2006]|Give'on et al.|2007l[Pueyo et al.|2009| ) are 

not suitable to find DM shapes that mitigate apertures 
discontinuities. We circumvent these numerical limita- 
tions by calculating DMs shapes that are based on the 
full non-linear problem, right panel of Fig. [T] 

3.3. The SR-Fresnel approximation 

In the general case, starting from Eq. |13| and following 
the derivation described in §.[5j the fieldat DM2 can be 
written as follows: 



E DM 2{x2,y 2 ) = |^dct^^^^ 1 (e,77)e ^2 " (^/l+Wl) e- ^ ^(^ ^2 + ^ , ' 2 ) d&fojj 



(16) 



02,2/2) 



where -Edmi (£,??) is the Fourier transform of the tele- 
scope aperture, TV and stands for the Fourier plane. 
We call this integral the Stretched- Remapped Fresnel ap- 
proximation (SR-Fresnel). Moreover det[J(x 2 , y 2 )] is the 
determinant of the Jacobian of the change of variables 
that maps (x2,y 2 ) to (xi,yi): 



det[J(x 2l y2)\ 



dfidgx 
dx dy 



dgi 

dx 



(17) 



02,2/2) 



r 



In the ray optics approximation, A ~ 0, the non linear 
transfer function between the two DMs becomes: 



(18) 

02,2/2) 



E D M2{x 2 ,y 2 ) = I \/\ det[J]\E DM i (/i,5i)} 
[E D M2{x 2 ,y2)f ! = |det[J] [E DM1 (/i,ffi)] 2 | 



(19) 



02,2/2) 



The square form (e.g Eq. 19 ) of this transfer function can 
also be derived based on conservation of energy princi- 
ples and is a generalization to arbitrary geometries of the 
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Figure 3. Schematic of the optical system considered: the telescope apertures is followed by two sequential Defromable Mirrors (DMs) 
in non-conjugate planes whose purpose is to remap the pupil discontinuities. The beam then enters a coronagraph to suppress the bulk of 
the starlight: in this figure we show an Apodized Lyot Pupil Coronagraph (APLC). This is the coronagraphic architecture we consider for 
the remainder of the paper but we stress that the method presented herein is applicable to any coronagraph. 



equation driving the d esign of PIAA coronagraphs ( Van- 
derbei & Traub||2005|). A full diffractive optimization of 



the DM surfaces requires use of the complete transfer 
function shown in Eq. |16| However, there do not exist 
yet tractable numericalmethod to evaluate Eq. [16] effi- 
ciently enough in order for this model to be included in 
an optimization algorithm. Moreover even solving the 
ray optics problem is extremely complicated: it requires 
to find the mapping function (fi,gi) which solv es the 
non-linear partial differential equat ion i n Eq. |19| Sub- 
stituting for(/i,3i) and using Eqs. |12b| yields a second 
order non-linear partial differential equation in This 
is the problem that we set ourselves to tackle in the next 
section, and is the cornerstone of our Adpative Com- 
pensation of Aperture Discontinuities. As a check, one 
can verify that in the small deformation regime (e.g. if 
h\{x, y) = Aecos(^j (mx+ny)) and h 2 (x,y) = —hi{x,y)) 



Eq. [19] yields the s ame phase-to -amplitude coupling as 
m Talbot imaging ( |Pueyo||2 008 [) . Eq.|19 | is a well know 
optimal transport problem (Mongc 17ST), which has al- 
ready been ide ntified as underlying op tical illumination 
optimizations ( |Glimm fc Ol ikcr 2002). While the exis- 
tence and uniqueness of solutions in arbitrary dimensions 
have been e xtensively discussed in the m athematical lit- 
erature (see Dacorogna & Moser ([1990) for a review), 
there was no practical numerical solution published up 
until recently. In particular, to our knowledge, not even 
a dimensional solution for which the DM surfaces can be 
described using a realistic basis-set has been published 
yet. We now introduce a method that calculates solu- 



tions to Eq. 19 which can be represented by feasible DM 
shapes. 

4. CALCULATION OF THE DEFORMABLE MIRROR 
SHAPES 

4.1. Statement of the problem 

Ideally, we seek DM shapes that fully cancel all the 
discontinuities at the surface of the primary mirror and 
yield a uniform amplitude distribution, as shown in the 
top panel of Fig. [4] A solutions for a particular geom- 
etry with four se condary support has been derived by 
Lozi et al. ( 2009[) . It relies on reducing the dimension- 
ality of the problem to the direction orthogonal to the 
spiders. It is implemented using a transmissive correct- 
ing plate that is a four-faced prism arranged such that 
the vertices coincides with the location of the spiders. 
The curvature discontinuities at the location of the spi- 
ders are responsible for the local remapping that removes 
the spiders in the coronagraph pupil. However such a so- 
lution cannot be readily generalized to the case of more 
complex apertures, where the secondary support struc- 
tures might vary in width, or in the presence of segment 
gaps. Moreover it is transmissive and thus highly chro- 
matic. Here we focus on a different class of solutions 
and seek to answer a different question. How well can 
we mitigate the effect of pupil discontinuities using DMs 
with smooth surfaces, a limited number of actuators (e.g 
a limited maximal curvature) , and a limi ted stroke? Un- 
der these constraints directly solving Eq. 19 (e.g. solving 
the forward problem illustrated in the top panel of Fig.[4J) 
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The forward problem 




The tapered reverse problem 




A 



Figure 4. Top: In the ideal case the two DMs would fully remap 
all the discontinuities in the telescope's aperture to feed a fully 
uniform beam to the coronagraph. However this would require 
discontinuities in the mirror's curvatures which cannot be achieved 
in practice. Moreover solving the Monge-Ampere Equation in t his 
direction is a difficult exercise as the right hand side of Eq, |19| 
presents an implicit dependence on the solution h'2- Bottom: We 
circumvent this problem by solving the reverse problem, where 
the the input beam is now uniform and the implicit dependence 
drops out. Moreover we taper the edges of the discontinuities by 
convolving the target field A(x\ ,y\) by a gaussian of full width at 
half maximum ui (ui = 50 cycles per aperture in this figure) . 

is not tractable as both factors on the left hand side of 
Eq. [19] depend on h 2 . More specifically, the implicit de- 
pendence of Ed mi (/i(^2,2/2),ffi(a;2,2/2)) on h 2 can only 
be addressed using finite elements solvers, whose solu- 
tions might not be realistically representable using a DM. 
However this can be circumvented using the reversibility 
of light and solving the reverse problem, where the two 
mirrors have been swapped. Indeed, since 



x-2 = h{h{x2,y2)>gi{x2,y2)) 



(20) 
(21) 



then we have the following relationship between the de- 
terminants of the forward and reverse remappings: 



1 



d/2 dg 2 
dx dy 



dg 2 
dx 



dfi dg 1 
dx dy 



dgi 

dx 



Moreover we are only interested in compensating asym- 
metric structures located between the secondary and the 
edge of the primary. We thus only seek to find (f 2 ,g 2 ) 
such that: 



E D Mi{xi,yi) = A(xi,yi) = 



[P(x,y)-{1-P (x,y))]*e 



(25) 



where Po{ x \tDi) is the obscured pupil, without seg- 
ments or secondary supports. Finally, we focus on so- 
lutions with a high contrast only up to a finite OWA. 
We artificially taper the discontinuities by convolv- 
ing the control term in the Monge Ampere Equation, 
[P(x,y) — (1 — Po{x,y))}, with a gaussian of width lu. 
Note that this tapering is only applied when calculating 
the DM shapes via solving the reverse problem. When 
the resulting solutions are propagated through Eq. [19] 
we use the true telescope pupil for Er>ui(xi, yi)- The 
parameter to has a significant impact on the final post- 
coronagraphic contrast. Indeed we are here working with 
a merit function that is based on a pupil-plane residual, 
while ideally our cost function should be based on image- 
plane intensity. By convolving the control term in the 
reverse Monge Ampere Equation, we low pass filter the 
discontinuities. This is equivalent to giving a stronger 
weight to low-to-mid spatial frequencies of interest in the 
context of exo-planets imaging . For each case presented 
in §.[6] we calculate our DM shapes over a grid of values 
of tu which correspond to low pass filters with cutoff fre- 
quencies ranging from ~ OWA to ~ 2 OWA and we keep 
the shapes which yield the best contrast. This somewhat 
ad-hoc approach can certainly be optimized for higher 
contrasts. However such an optimization is beyond the 
scope of the present manuscript. 

The problem we are seeking to solve is illustrated in 
the second panel of Fig. [4] In this configuration the full 
second order Monge Ampere Equation can be written as: 



(1- 



z d2hl ) 



dxdy J 



A(x,y) 2 (26) 



where we have dropped the (xi,y±) dependence for clar- 
ity. Since we are interested in surface deformations which 
can realistically be created using a DM, we seek for a 
Fourier representation of the DMs surface: 



02,1/2) 



(22) 

We thus focus on the inverse problem, bottom panel of 
Fig. H that consists of first finding the surface of hi as 
the solution of: 



h 1 (x,y)= I ^H 1 (X,Y) 



~Z 



N/2 N/2 
n=-N/2 m=-N/2 



i^g-(mX+nX) 



(27) 



<^[E DM 2{f2,g2) 

= [E D Mi(xi,yi)} 2 



dh dg 2 
dx dy 



dg 2 
dx 



with a_ 



a 



where N is the limited num- 



(23) 

Since our goal is to obtain a pupil as uniform a possible 
we seek a field at DM2 as uniform as possible: 



ber of actuators across the DM. Note that we have nor- 
malized the dimensions in the pupil plane X = x/D, 
Y = y/D. The normalized second order Monge Ampere 
Equation is then: 



E D M2(f2(xi,yi),92(xi,y 1 )) = 

J J Eomi(x, y) 2 dxdy — Constant. 



(1 



d 2 Hi 
dX 2 



)(1 



d 2 #i. 
dY 2 ' 



( d 2 H, V 
l dXdY ) 



A(X,Y) 2 (28) 



(24) 



For each configuration in this paper we first solve Eq. [28] 
and then transform the normalized solution in physical 
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Forward coordinate system 

Hi(X,Y) 




Reverse coordinate system 



HiiX.Y) 




H 2 (X,Y) 




Figure 5. Forward and reverse coordinate systems, with their 
respective forward and reverse coordinate transforms, in the case 
of the folded system studied in this paper. Note that this figure 
shows normalized units. As explained in the body of the text the 
correspondence between normalized and real units scales as follow: 

favi) = [DXi,DYi), {fi,9i) = (DFi, DGi), H t = ^-h h where 
i = 1,2, D is the aperture diameter and Z the separation between 
DMs. We solve Monge Ampere Equation in normalized coordinates 
and then apply the scalings in order to find the true DM shapes. 



units, which depends on the DMs diameter D and their 
separation Z . 

4.2. Solving the Monge- Ampere equation to find Hi 

Over the past few years a nu mber of numerical algo- 
rithms a imed at solving Eq. 28 h ave emerged in the lit 



erature (Loeper & Rapetti 2005 Benamou et al. 



2010] ) . 

Here we summarize our imple mentation ot two of them : 
an explicit Newton method jLoeper fe Rapetti 20051 ), 
and a semi-implicit method ( |Froese fe Uberm an 2012). 
We do not delve into the pr oof of convergence of eac h 
method, they can be f o und in|Loeper fe R apetti (2005); 

ote 



Benamou et al.| (|20lot: iFroese fc Uberman| ( 20T2| 
that Zheligovsky et al. Q2010P discussed both approaches 
in a cosmological context and devised Fourier based solu- 
tions. Here we are interested in a two dimensional prob- 
lem and we outline below the essence of each algorithm. 



4.2.1. Explicit Newton algorithm 

This method was first introduced by|Loeper fc Rapetti| 
( 2005 1 and relies on the fact that Eq. 28 can be re- written 



as 



det 



1 + 



, dX 2 
d 2 H t 
dXdY 



d 2 H ± 
dXdY 

1 + dY 2 



= det[Id+H(H 1 (X, Y))} 
(29) 

where %(•) is the two dimensional Hessian of a scalar field 
and Id the identity matrix. If one writes H\ — u + v with 
IMI <C Hull then: 



J 



det[Id+rl{u + Sv)} = det [Id + H(u)]+ S Tr {Id + U{u)) ] H(v) + o{5 2 ) 



where (-)t denotes the transpose of the comatrix. Eq. 28 

d 2 u . d 2 v „ d 2 u 



(1 



d 2 u d 2 v 

Wy^'ox 2 



(1 



dX 2 ' dY 2 dXdY dXdY 



can thus be linearized as: 

' A{X,Y) 2 ~det[H( 



X 2 + Y 2 



(30) 



(31) 



The explicit Newton algorithm relies on Eq. 31 and can then be summarized as carrying out the following iterations: 

• Choose a first guess H®. 

• At each iteration k we seek for a solution of the form H± +1 = H\ + V k , where V k is the DM shape update. 

• In order to find V k we write: 



(1 



d 2 H k d 2 V k ^d 2 H k d 2 V k 



dX 2 ' dY 2 



dXdY dXdY 



Re<M) = -(A 2 - det[Id + U{H k )]) 



(32) 
(33) 



and solve 



Eq. 



38 



is a linear partial differential equation in V 
in a Fourier series we write V k as: 

N/2 

E 

= -N/2i 



L E (H k ,V k )=R E (H k ). (34) 
/; Since we are interested in a solution which can be expanded 



V k (X,Y) 



N/2 

E 

=-N/2 



v k ( 



i^imX+nX) 



(35) 



Both the right hand side and the left hand side of Eq. [38] can be written as a Fourier series, with a spatial 
frequency content between — N and N cycles per aperture. Equating each Fourier coefficient in these two series 
yields the following linear system of (2N + l) 2 equations with (N + l) 2 unknowns. 
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For all too, no £ {—N,N]: [ e ^(m x+n Y) ^ g ^ yk^ _ R E (H k )] dX dY = (36) 

J DM 

When searching for V k as a Fourier series over the square geometry chosen here, this inverse problem is always 
well posed. 

• Update the solution H% +1 = + V k+1 



The convergence of this algorithm reli es on the introduc- 



tion o f a damping constant r > 1. Loeper & Rapetti 

(2005) showed that as long as x + r + H k remains 
convex, which is always true for ACAD with reason- 
ably small aperture discontinuities, there exists a r large 
enough so that this algorithm converge towards a solu- 
tion of Eq. 28 However since this algorithm is gradient 
based, it is not guaranteed that it converges to the global 
minimum of the underlying non-linear problem. In order 
to avoid having this solver stall in a local minimum we 
follow the methodology outlined by|Froese fc Oberman 



to get within a reasonable neighborhood of the global 
minimum. 

4.2.2. Implicit algorithm 

This algorithm, along with its convergence pro of, is 
thoroughly explained in Froese & Oberman ( 2012 ) . It 
relies on rewriting Eq. |28| as: 



(2012) and first carry out a series of implicit iterations 



d 2 H x (PHi 
dX 2 + dY 2 



= v/dct[/d + H(H 1 (X, Y))} 2 + 2A(X, Y) 2 

(37) 



The implicit method consists of carrying out the following iterations: 

• Choose a first guess 

• In order to find H k+1 we write: 

d 2 H^ 



Li (Hi, V k ) 



k+l 



dX 2 



d 2 H k 
dY 2 



R E (H k ) = v/dct[/d + %(Hi(X, Y))] 2 + 2A(X, Y f 



and solve 



Lj(H k+1 ) = Rj(H k ] 



(38) 



This problem is a linear system of (N- 
on a Fourier Basis: 

For all too, n € [-N,N]: 



T) equations with (N + l) 2 unknowns and can be solved using projections 



a i2ir(m X+n Y) 



DM 



[L / (i7 1 fc , V k ) - Ri{H k )} dX dY = 



(39) 



Note that the term under the square root in Rj(H k ) is guaranteed to be positive at each iteration. 
• Iterate over k 



The inverse problem in Eq. 39 is always well posed, for 
any basis set or pupil geometry, while the explicit New- 
ton method runs into convergence issues when not us- 
ing a Fourier basis over a square. When seeking to use 
a basis set that is more adapted to the geometry of the 
spiders and segments or when using a trial influence func- 
tion basis for the DM, the implicit method is the most 
promising method. In this paper we have limited our 
scope to solving the reverse problem in the bottom panel 
of Fig. [4 and used a Fourier representation for the DM, 
we are able to use both methods. In order to make sure 
that the algorithm converges towards the true solution of 
Eq. 28 we first run a few tens of iterations of the implicit 
method and, once it has converged, we seek for a more 
accurate solutions using the Newton algorithm. Typical 
results are shown in Fig. [6] where most of the residual 
error resides in the high spatial frequency content (e.g. 



I 

above N cycles per aperture). Our solutions are limited 
by the non- optimality of the Fourier basis to describe 
the mostly radial and azimuthal structures present in 
telescopes's apertures. Moreover the DM shape is the 
result of the minimization of a least squares residual in 
the virtual end-plane of the reverse problem, with little 
regard to the spatial frequency content of the solution 
in the final image plane of the coronagraph. While this 
method yields significant contrast improvements, as re- 
ported in §.[6j we discuss in §.[7] how it can be refined for 
higher contrast. 

4.3. Deformation of the Second Mirror 

Once the surface of DM1 has been calculated as a so- 
lution o f Eq. fl9| we compute the surface of DM2 based 



on Eqs. 12b which stem from enforcing flatness of the 
outgoing on-axis wavefront. We seek a Fourier represen- 



12 Pueyo & Norman 




Figure 6. Virtual field at Ml when solving the reverse problem. A(Xi, Yi) is the desired apodization (top), A n (X±, Y\) is the apodization 
obtained after solving the Monge Ampere Equation (center). The bottom panel shows the difference between the two quantities: the bulk 
of the energy in the residual is located in high spatial frequencies that cannot be controlled by the DMs. 





' " (\-F 2 ( 1 -,Y 1 ) 



X-2 — Fi (X i , Y\ ) 




„ 

$ - 



Figure 7. Boundary conditions seen in the horizontal remapped 
space. Should the boundary conditions have strictly been enforced 
by our solver then F 2 (-\,Y) = -1, F 2 {\,Y) = §, F 2 (X,-\) = 
X, F 2 {X, §) = X. The remapping function obtained with our 
solutions are very close to these theoretical boundary conditions 
and the residuals can easily be mitigated by sacrificing the edge 
rows and columns of actuators on each DM. 



tation for the surface of DM2: 

h 2 (x,y) = ^-H 2 (X,Y) = 

N/2 N/2 

^ E E 6 m ,„e^(™™> (40) 

n=-N/2 m=-N/2 

Plugging the s olution found in the previous step for h\ 
into Eqs. |llb| yields a closed form for the normalized 
remapping functions, (F 2 ,G 2 ): 

N/2 N/2 

F 2 (X ll Y 1 ) = X 1 - ]T ]T flma™^ 12 '^^' 

n=-N/2 m=-N/2 
N/2 N/2 

G 2 (X 1 ,Y 1 )=Y 1 - ]T a™ a m , n e l2 ^ mX i +nY ^ 

n=-N/2 m=-N/2 

Then the normalized version of Eqs. |12b"| can be rewritten 
as: 

N/2 

L^XuYt) = i2lTm e^F^Y^+nF^Y,)) 

m=-N/2 

R X (X 1 ,Y 1 )=X 1 -F 2 (X 1 ,Y 1 ) 

L X = R X (41) 

N/2 
m=-N/2 

R y (X 1 ,Y 1 )=Y 1 -G 2 (X 1 ,Y 1 ) 

L y = Ry (42) 
We then multiply each side of Eq. [4l] and Eq. [42] by: 

^(moF.iXuY^+noF^XuY,)) det [ Id + U (X^))] 

where (mo, no) corresponds to a given DM spatial fre- 
quency. 



Integrating over the square area of the DM and using the orthogonality of the Fourier basis yields the following 
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system of 2 * (N + l) 2 equations with (N + l) 2 real unknowns: 
For all (too, no): 

2m m b mQino = f R X (X, Y) det[7d + H(H^(X, Y))]e l2 * {maF * {x ' Y)+n ° F ^ x > Y)) dXdY 

J DM 

For all (mo, tlq): 

Vim n b mo . no = f R X {X, Y) det[/d + H(H?(X, Y))]^ 2 < m ° F ^ x ^+ n ° F ^ x ' Y ^ dXdY 

J DM 

We then find H^ , the normalized surface of DM2, by solving this system in the least squares sense. 

I 




Figure 8. Boundary conditions seen in the vertical remapped 
space. Should the boundary conditions have strictly been enforced 
by our solver then G 2 (-\,Y) = Y, G 2 {\,Y) = Y, G 2 (X,-A) = 
— G 2 (X, |) = \- The remapping function obtained with our 
solutions are very close to these theoretical boundary conditions 
and the residuals can easily be mitigated by sacrificing the edge 
rows and columns of actuators on each DM. 

Once the Monge Ampere Equation has been solved, the 
calculation of the surface of the second mirror is a much 
easier problem. Indeed, by virtue of the conservation 
of the on-axis optical path length, finding the sur face of 
DM2 only consists of solving a linear system (see Traub 
fc Vanderbei|([2003|). 



4.4. Boundary Conditions 

The method described above does not enforce any 
boundary conditions associated with Eq. |28| One prac- 
tical set of boundary conditions consists of forcing the 
edges of each DM to map to each other: 



F l (X,±\)=X 
G i (±\,Y)=Y 



(43) 
(44) 
(45) 
(46) 



with i = 1,2. These correspond to a set of Neumann 
boundary conditions in Hi (X, Y) and H2 (X, Y) . These 
boundary conditions can be enforced by a ugm enting the 
dimensionality of the linear systems on Eq. |36| and Eq. |39[ 
however doing so increases the residual least squares er- 
rors and thus hampers the contrast of the final solution. 
Moreover Fig. [7] and Fig. [8] show that, because of the one 
to one remapping near the DM edges in the control term 
of the reverse problem, the boundary conditions are al- 
most met in practice. For the remainder of this paper we 
thus do not include boundary conditions when calculat- 



ing the DM shapes, when solving for Hi(X, Y) in Eq. 128 
since, in the worse case, only the edge rows and columns 
of the DMs actuator will have to be sacrificed in order 
for the edges to truly map to each other. 

4.5. Remapped aperture 

For a given pupil geometry we have calculated 
(Hi, )■ We then convert the DM surfaces to real units, 
(hi,h2), by multiplication with D 2 / Z. We e valuate the 
remapping functions using Eqs. | lib] and [l2b| and obtain 
the field at the entrance of the coronagraph in the ray 
optic approximation 



E 



DM2 



(X2,V2) - {v^M^Mi(/i,9i)e i "^ 1 ^)+' 11 ^^)-' 12 )} 



(47) 



where the exponential factor corresponds to the 
Optical Path Length through the two DMs. Even if 
the surface of the DMs has been calculated using only 
a finite set of Fourier modes, we check that the optical 



r 



path length is conserved. 



Fig 
limitcc 



shows that since the curvature of the DMs is 
jy the number of modes N, our solution docs 
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Before DMs After DMs 





Figure 9. Beam amplitude before and after the DMs in the case 
of a geometry similar to JWST. We chose to solve the reverse 
problem over a square using a Fourier basis-set and not stricly 
enforcing boundary conditions. This results in small distortions of 
the edges of the actual aperture in the vicinity of segments gaps and 
secondary supports. We address this problem by slightly oversizing 
the secondary obscuration and undersizing the aperture edge in the 
coronagraph. 



not fully map out the discontinuities induced by the sec- 
ondary supports and the segments. However, they are 
significantly thinner and one can expect that their im- 
pact on contrast will be attenuated by orders of magni- 
tude. In order to quantify the final coronagraphic con- 
trasts of our solution we then propagate it through an 
APLC coronagraph designed using the method in §. [2j 
In the case of a hexagon based primary (such as JWST), 
we use a coronographic apodizer with a slightly over- 
sized secondary obscuration and undersize outer edge in 
order to circularize the pupil. Note that this choice is 
mainly driven by the type of coronagraph we chose in 
§. [2] to illustrate our technique. Since the DM control 
strategy presented in this section is independent of the 
coronagraph, it can be generalized to any of the starlight 
suppression systems which have been discussed in the 
literature. For succinctness we present our results using 
coronagraph solely based on using pupil apodization (ei- 
ther in an APLC or in a PIAAC configuration). Results 
for a JWST geometry are shown on Fig. [11] and Fig. [12 



and discussed in §. [6] Eq. [47| assumes that the propaga- 
tion between the two DMs occurs according to the laws 
of ray optics. In the next sect ion we derive the actual 
diffractive field at DM2, e.g Eq. [l6j and show that in the 
pupil remapping regime of ACAD, edge ringing due to 
the free space propagation is actually smaller than in the 
Fresnel regime. 

5. CHROMATIC PROPERTIES 

5.1. Analytical expression of the diffracted field 

ACAD is based on ray optics. It is an inherently broad- 
band technique, and provided that the coronagraph is 
optimized for broadband performance ACAD will pro- 
vide high contrast over large spectral windows. However, 
when taking into account the edge diffraction effects that 
are captured by the quadratic integral in Eq. |16[ the 
true propagated field at DM2 becomes wavelength de- 
pendent. More specifically, when A is not zero then the 
oscillatory integral superposes on the ray optics field a 
series of high spatial frequency oscillations. In theory, 
it would be best to use this as the full transfer func- 
tion to include chromatic effects in the computation of 
the DMs shapes. However, as discussed in S. HJ solv- 
ing the non-linear Monge- Ampere Equation is already a 
delicate exercise, and we thus have limited the scope of 
this paper to ray optics solutions. Nonetheless, once the 
DMs' shapes have been determined using ray optics, one 
should check whether or not the oscillations due to edge 
diffraction will hamper the contrast. This approach is 
reminiscent of the design of PIAA systems where the mir- 
ror shapes are calculated first using geometric optics and 
are then propagated through the diffractive integral in 
order to check a posteriori whether or not the ch romatic 
diffractive artif acts are below the design contrast ( Pluzh- 
nik et aL| "2005). In this section we detail the derivation 
of Eq. |l5[that is the diffractive integral for the two DMs 
remapping system and use this formulation to discuss the 
diffractive properties of ACAD. 



We start with the expression of the second order diffractive field at DM2 as derived in Pueyo et al. (2011), Eq. 13 



Edm2(x2,V2) = < / E DM i(x,y)e' 



^ {x - Xl f+2 l ^( X ~x 1 )(y~y 1 )+ 3 l ^(y-y 1 f 



dxdy 



We write Edmi{ x i y) as its inverse Fourier transform 

E DM i(x,y) = [ E DM1 {i,r,)e i2 <^ + y^ d£dr, 



(48) 



and insert this expression in Eq. |13| Completing the squares in the quadratic exponential factor then yields: 



E D M2{X2,D2) = 



iXZ 



^Z5Mi(^^)/i(/i,.9i)e^ [/l « +9l ^ 1 e^ AZ (^ | (^'«i'« 2+ ^ l <^-")' )2 ) dftfy 



(49) 



>2,2/2) 



with 



^ {x - Xl -^f+2^{ X ~ Xl )(y-y 1 )+°^(y-y 1 -^f 



dxdy 



(50) 



(zi>?/i) 



The integral over space, Ii(y\,x{), is the integral of a complex gaussian and can be evaluated analytically. This 
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yields: 



E D M2{X2,V2) 



I §_h dg 2 

dx dy 



^Mi(^)e ,2lIlie+Ml e-" AZ (^^^'' 2 )^ 



(51) 



We thus have expressed -EpA^fei y 2 ) as a function of [x\, yi) — (/i(x2j2/2)j5i(2 ; 2jZ/2))- This expression can be further 

Hone can derive ^£\ {xim) = ^ IT L ^ ^\(,,,^ and ^ - - 1 ^ 

rel 



simplified: using Eqs. 
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Which finishes to prove 



to22 

Eq, 



J TV 



dct[J(x 2 ,y 2 )} dy l(>2, yi) dncl dy K*i>»l) — det[J(z2,i/2)] 9a: 1(^2, 2/2)- 

(52) 



(X2,V2) 



This expression is very similar to a modified Fresnel propagation and can be rewritten as such: 

E DM2 (x 2 ,y 2 ) = {^ DM1 (x, 2 /)e-^(^)- 1 ^^) 2+ ^)- 2 ^) 2 ) dxdy 



(53) 



(X2,y2) 



Because of this similarity we call this integral the Stretched- Remapped Fresnel approximation (SR- Fresnel). Indeed 

s 91 dfi 

in this approximation the propagation distance is stretched by ( g^fji , g^Jjj ) and the integral is centered around the 
remapped pupil (A, 51). 

r 



5.2. Discussion 

The integral form provides physical insight about the 
behavior of the chromatic edge oscillations. If we write 



T x = 



r, 



det[J] 



dgi 

dy 

det[J] 



(x2,y 2 ) 



dx 



(54) 



(55) 



(3:2,1/2) 

we can identify a several diffractive regimes: 

1. when T x — T y = 1 the £1)4/2(12, 2/2) reduces to a 
simple Fresnel propagation. Edge ringing can them 
be mitigated using classical techniques such as pre- 
apodization, or re-imaging into a conjugate plane 
using oversized relay optics. 

2. when T x ,T y < 1 then it is as if the effective propa- 
gation length through the remapping unit was in- 
creased. This magnifies the edge chromatic ringing. 

3. when T xl T y > 1 then it is as if the effective propa- 
gation length through the remapping unit was de- 
creased. This damps the edge chromatic ringing. 

4. when T x > 1 and T y < 1, e.g. at a saddle point in 
the DM surface, then it is as if the effective propa- 
gation length through the remapping unit was de- 
creased in one direction and increased in the other. 
The edge chromatic ringing can either be damped 
or magnified depending on the relative magnitude 
of T x and T y . 

In the case of a PI A A coronagraphs, the mirror shapes 
are such that T x ,T y > 1 at the center of the pupil and 
, T y <C 1 at the edges of the pupil, where the discon- 
tinuities occur. As a consequence the edge oscillations 
are largely magnified when c omp ared to Fresnel oscilla- 
tions (see right panel of Fig. 10 1, and apodizing screens 



are necessary in order to reduce the lo cal curvature of 
the m i rror's shape ( as als o discussed in Pluzhnik et al. 
( |2005| ); |Pueyo et al-lpm] )). In the case of ACAD, where 
the x-axis is chosen to be perpendicular to the disconti- 
nuity, the surface curvature is such that T x > 1, T y ~ 1 
at the discontinuities inside the pupil and F x > 1, T y ~ 1 
elsewhere. This yields damped chromatic oscillations at 
the remapped discontinuities and Fresnel osci llat ions at 
the edges of the pupil (see right panel of Fig. 10 1. Note 
that Fig. |10| was generated using a one dimensional toy 
model that assumes Eq. [TBI is s eparable, e.g T y = 1, as 



Pueyo et al.p|2t)ll[ ). In practice at the sad- 



described in 

die points of the optical surfaces, near the junction of 
two spiders for instance, j x > 1, T y <^ 1 and thus 
our separable model does not guarantee than in the true 
2D case chromatic edge oscillations might not be locally 
amplified. However even near the saddle points ACAD 
provides a strong converging remapping in the direction 
perpendicular to the discontinuity and very little diverg- 
ing re-mapping in the other direction. As a consequence 
r^, ^ 1 and r y is smaller than but close to one. We thus 
predict that even at these locations chromatic ringing will 
not be amplified. Even if ACAD based on pupil remap- 
ping, its diffraction properties are qualitatively very dif- 
ferent from PIAA coronagraphs since edge ringing is not 
amplified beyond the Fresnel regime at the pupil edges, 
and is attenuated near the discontinuities. We conclude 
that in most cases ACAD operates in a regime where edge 
chromatic oscillations are not larger than classical Fres- 
nel oscillations, and sometimes actually smaller. As a 
consequence the chromaticity of this ringing can be miti- 
gated using standard techniques developed in the Fresnel 
regime and we do not expect this phenomenon to be a 
major obstacle to ACAD broadband operations. 

5.3. Diffraction artifacts in ACAD are no worse than 
Fresnel ringing 

We have established that the diffractive chromatic os- 
cillations introduced by the fact that DM1 and DM2 are 
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Figure 10. Left: Comparison of edge diffraction between PIAA and Fresnel propagation. Because the discontinuities in the pupil 
occur at the location where F x < 1, a PIAA amplifies the chromatic ringing when compared to a more classical propagation. Right: 
Comparison of edge diffraction between ACAD and Fresnel propagation. Because the discontinuities in the pupil occur at the location 
where > 1, a ACAD damps the chromatic ringing when compared to a more classical propagation. These simulations were carried out 
for a one dimensional aperture. In this case the Monge Ampere Equation can be solved using finite elements and the calculation of the 
diffractive effects reduces to the evaluation of various Fresnel special functions at varying wavelength. The full two-dimensional problem is 
not separable and requires the development of novel numerical tools. 

• The coronagrahic apodizer can be placed in a plane 
that is conjugate to the DM1. This can be achieved 
by re-imaging DM2 through a system of oversized 
optics (the over-sizing factor increases steeply when 
the pupil diameter decreases). By definition there 
are no Fresnel edge oscillations in such a plane. Al- 
ternatively a coronagraph without a pupil apodiza- 
tion (amplitude or phase mask in the image plane) 
can be used, and in this configuration it is only 
sufficient for the optics to be oversized. 



not located in conjugate planes is no worse than classi- 
cal Fresnel ringing from the aperture edges and can be 
mitigated using well-know techniques which have been 
developed for this regime. While a quantitive tradeoff 
study of how to design a high contrast instrument which 
minimizes such oscillations regime is beyond the scope of 
this paper, we briefly remind their qualitative essence to 
the reader: 

• The edges of the discontinuities in the telescope 
aperture can be smoothed via pupil apodization 
before DM1. This solution is not particularly ap- 
pealing as it requires the introduction of a trans- 
missive, and thus dispersive, component in the op- 
tical train. 

• The distance between the two DMs can be reduced. 
Indeed the DMs deformations presented herein, for 
3 cm DMs separated by 1 m, are all < 1 /im while 
current technologies allow for deformations of sev- 
eral microns. As the edge ringing scales as Z/D 2 
chromatic oscillations will be reduced by decreasing 
Z. Since the DM surfaces scale as D 2 /Z reducing 
Z will increase the DM deformations but have little 
impact on the feasibility of our solution as current 
DM technologies can reach 4 /im strokes. 



Note that these three solutions are not mutually exclu- 
sive and that only a full diffractive analysis, which uses 
robust numer ical propagators that have been developed 
based on Eq. |16[ can quantitatively address the trade- 
offs discussed above. The deve lopment of such pro pa- 



gators is our next priority. In Pueyo et al. (2011) we 



laid out the theoretical foundations of such a numerical 
tool in the case of circularly symmetric pupil remapping 
and this solution has be en since then pract ically imple- 
mented, as reported by Krist et al. (20101. Generaliz- 
ing this method to a tractable propagator in the case 
of arbitrary remapping is a yet unsolved computational 
problem. In the meantime we emphasize that while the 
spectral bandwidth of coronagraphs whose incoming am- 
plitude has been corrected using ACAD will certainly be 
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Surface of DM1 Surface of DM2 




Figure 11. Results obtained when applying our approach to a geometry similar to JWST . We used two 3 cm DMs of 64 actuators 
separated by 1 m. Their maximal surface deformation is 1.1 fim, well within the stroke limit of current DM technologies. The residual 
light in the corrected PSF follows the secondary support structures and can potentially be further cancelled by controlling the DMs using 
an image plane based cost function, see Fig. |23| 
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Correction with N - 64 

DMs at rest 

■ No spiders 




A/D in case of an APLC 

Figure 12. Case of JWST: Radial average obtained when ap- 
plying ACAD. We used two 3 cm DMs of 64 actuators separated 
by 1 m. Their maximal surface deformation is 1.1 fim, well within 
the stroke limit of current DM technologies. ACAD yields a gain 
in contrast of two orders of magnitude, and provides contrasts lev- 
els similar to upcoming Ex-AO instruments, which are designed 
on much friendlier apertures geometries. Since ACAD removes the 
bulk of the light diffracted by the asymmetric aperture disconti- 
nuities, the final contrast can be improved by contr ollin g the DMs 
using and image plane based cost function, see Fig. |23| 



limited by edge diffraction effects, but these effects are no 
worse than Fresnel ringing and can thus be mitigated us- 
ing optical designs which are now routinely use d in high 
contrast instruments (see Verinaud et al. (2010) for such 
discussions) . For the remainder of this paper we thus as- 
sume the diffractive artifacts have been adequately miti- 
gated and we compute our results assuming a geometric 
propagation between DM1 and DM2. 



6. RESULTS 
6.1. Application to future observatories 
6.1.1. JWST 

We have illustrated each step of the calculation of the 
DM shapes in §. fusing a geometry similar of JWST. 
This configuration is somewhat a conservative illustra- 
tion of an on-axis segmented telescope as it features thick 
secondary supports and a "small" number of segments 
whose gaps diffract light in regions of the image plane 
close to the optical axis (the first diffraction order of a 
six hexagons structure is located at ~ 3X/D). In order 
to assess the performances of ACAD on such an obser- 
vatory architecture we chose to use a coronagraph de- 
signed around a slightly oversized secondary obscuration 
of diameter 0.25 D, with a focal plane mask of diameter 
8 A/D, an IWA of 5 X/D and an OWA of 30 A/D. The field 
at the entrance of the coronagraph after rema ppin g by 
the DMs is shown on the top right panel of Fig. [IT] The 
DM surfaces, calculated assuming 64 actuators across the 
pupil (N — 64 in the Fourier expansion) and DMs of di- 
ameter 3 cm separa ted by Z = 1 m, are shown on the 
middle panel of Fig. 11 They are well within the stroke 
limit of current DM technologies. The surfaces were cal- 
culated by solving the reverse problem over an even grid 
of 10 cutoff low-pass spatial frequencies ranging between 
30 and 70 cycles per apertures for the tapering kernel oj. 
The value yielding the best contrast was chosen. Note 
that the optimal cutoff frequency depends on the spa- 



tial scale of the discontinuities, and that higher contrasts 
could be obtained by choosing a set of two convolution 
kernels for the reverse problem and finding the optimal 
solution using a fin er g rid. However, the results in the 
bottom row of Fig. [IT] are extremely promising. Fig. [12] 
shows a contrast improvement of a factor of 100 when 
compared to the raw coronagraphic PSF, which is quite 
remarkable for an algorithm which is not based on an 
image-plane metric. These results illustrate that even 
with a very unfriendly aperture similar to JWST one 
can obtain contrasts as high as envisioned for upcoming 
Ex-AO instruments, which have been designed for much 
friendlier apertures. While we certainly do not advocate 
to use such a technique on JWST, this demonstrates that 
ACAD is a powerful tool for coronagraphy with on-axis 
segmented apertures. 

6.1.2. Extremely Large Telescopes 

We now discuss the case of Extremely Large Telescopes 
and provide an illustration using the example of the 
Thirty Meter Telescope. We considered the aperture ge- 
ometry shown on the top left panel of Fig. [13] It consists 
of a pupil 37 segments across in the longest direction and 
a secondary of diameter ~ 0.1 2D which is held by three 
main thick struts and six thin cables. As seen on the 



bottom left panel of Fig. 13 the impact of segment gaps 
is minor as they diffract light beyond the OWA of the 
coronagraph. When using a coronagraph with a larger 
OWA the segment gaps will have to be taken into ac- 
count, and will have to be mitigated using DMs with a 
larger number of actuators. In order to obtain first order 
estimates of the performances of ACAD on the aperture 
geometry shown on the top left panel of Fig.[l3 we chose 
to use a coronagraph designed around a slightly over- 
sized secondary obscuration of diameter 0.15 D, with a 
focal plane mask of 6X/D diameter, an IWA of AX/D 
and an OWA of 30A/Z?. The field at the entrance of the 
coronagraph after remapp ing by the DMs is shown on 
the top right panel of Fig. [13] The DM surfaces, calcu- 
lated assuming 64 actuators across the pupil [N = 64 in 
the Fourier expansion) and DMs of diameter 3 cm sep- 
arated by Z = 1 m, are shown on the middle panel of 
Fig. [13] They are well within the stroke limit of current 
DM technologies. The surfaces were calculated by solv- 
ing the reverse problem over an even grid of 10 cutoff 
low-pass spatial frequencies ranging between 30 and 70 
cycles per apertures for the tapering kernel u). The value 
yielding the best contrast was chosen. The final PSF 
is shown on the bottom right panel of Fig. [13] and fea- 
tures a high contrast dark hole with residuaTdiffracted 
light at the location of the spiders' diffraction structures. 
The impact on coronagraphic co ntrast of secondary sup - 



ports was thoroughly studied by|Martinez et al. (2008 



They concluded that under a 90% Strehl ratio, the con- 
trast in most types of coronagraphs is driven by the sec- 
ondary support structures to levels ranging from 10~ 4 to 
10~ 5 . This, in turn, leads to a final contrast after post- 
processing (called Differential Imaging) of ~ 10~ 7 — 10~ 8 . 
Fig. [14] shows that using ACAD on an ELT pupil yields 
contrasts before any post-pr ocessing which are com para- 



ble to the ones obtained by Martinez et al. (2008) after 



Differential Imaging. This demonstrates that should two 
sequential DMs be integrated into a future planet finding 
instrument, setting their surface deformation according 
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Figure 13. Results obtained when applying our approach to a TMT geometry. We used two 3 cm DMs of 64 actuators separated by 1 
m. Their maximal surface deformation is 0.9 /im, well within the stroke limit of current DM technologies. The final contrast is below 10 7 , 
in a regime favorable for direct imaging of exo-planets with ELTs. 
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Figure 14. Case of TMT: radial average obtained when apply- 
ing ACAD. We used two 3 cm DMs of 64 actuators separated by 
1 m. Their maximal surface deformation is 0.9 /im, well within 
the stroke limit of current DM technologies. The final contrast is 
below 10 7 , in a regime favorable for direct imaging of exo-planets 
with ELTs. Since ACAD removes the bulk of the light diffracted 
by the asymmetric aperture discontinuities, the final contrast can 
be enhanced by controlling the DMs using and image plane based 
metric. 

to the methodology presenting above would allow this in- 
strument to perform its scientific program at a very high 
contrast. Moreover the surface of the DMs could be ad- 
justed to mitigate for the effect of missing segments at the 
surface of the primary (when for instance the telescope 
is operating while some segments are being serviced). 

6.2. Hypothetical cases 

6.2.1. Constant area covered by the secondary support 
structures 

In the case of ELTs with large number of small 
segments (when compared to the aperture size) , gaps 
diffract light far from the optical axis (see Fig. 13 for 
an example) . The secondary support structures are then 
the major source of unfriendly coronagraphic diffracted 
light. Under the assumption that thick structures are 
necessary to support the heavy secondary over the very 
large ELT pupils, one can use the aperture area covered 
by the spiders as a proxy of the secondary lift constraint. 
We have thus explored a series of geometries for which 
the number of spiders increases as they get thinner while 
the overall area covered by the secondary support struc- 
tures remains c ons tant. In the examples shown from 
Fig. [15] to Fig. [THJ the area covered by the secondary 
support structures is 1.5 times greater than in the TMT 
geometry discussed above. In all cases we used a corona- 
graph with a central obscuration of 0.15 D, with a focal 
plane mask of 6X/D diameter, and IWA of AX/D and an 
OWA of 30 A/D. The surfaces were calculated by solv- 
ing the reverse problem over an even grid of 10 cutoff 
low-pass spatial frequencies ranging between 30 and 70 
cycles per apertures for the tapering kernel. The value 
yielding the best contrast was chosen. This exercise leads 
to several conclusions pertaining to the performances of 
ACAD with various potential ELT geometries. 

Clocking of the spiders with respect to the DM 



of the clocking of the spiders with respect to the DMs ac- 
tuator grid (or the Fourier grid in our case). When the 
secondary support structures are clocked by 45° with re- 
spect to the DM actuators they are much more attenu- 
ated by ACAD, thus yielding higher contrast. This is 
an artifact of the Fourier basis set chosen and would 
be mitigated by using DMs whose actuator placement 
presents cir cular and azimuthal symmetries (Watanabe 
et al. 1120081). 



Annulus in the PSF with a large number of spi- 
ders 

When the number of secondary support struts becomes 
very large (> 20), an interesting phenomena occurs in the 
raw PSF: the spiders diffract light outside an annulus 
of radius Nspiders/n/D, just as spider web masks do in 
the ca se of shaped pupil coronagraphs (Vanderbei et al. 
2003b|. The "bump" located beyond that spatial fre- 



The top two panels of Fig. 15 illustrate the importance 



quency is more difficult to attenuate using the DMs (see 
Fig. 15 for an illustration). ACAD creates small ripples 
at the edges of the remapped discontinuities and when 
too many discontinuities are in the vicinity of each other, 
then these ripples interfere constructively and hamper 
the starlight extinction level yielded by ACAD. 

A lot of thin spiders is more favorable than a 
few thick spiders 

In general decreasing the width of the spiders while 
increasing their number is beneficial to the contrast ob- 
tained after ACAD as illustrated on the radial averages 
on Fig. 16 and Fig. 18 When one increases the number of 
spiders while decreasing their width in a classical corona- 
graph, the peak intensity of the diffraction pattern of one 
spider decreases as the squared width of the spider. The 
radially averaged contrast improvement without ACAD 
is then somewhat lesser than the square of the spider 
thinning factor as it is mitigated by the increasing num- 
ber of spiders. When using ACAD the spiders are seen 
by the coronagraph as much thinner than they actually 
are (by a factor r) and thus the peak intensity of their 
diffraction pattern is lower by a factor of r 2 . Our numer- 
ical experiments show that r increases when the spider 
width decreases. As consequence, the overall contrast 
gain after ACAD when decreasing the width of the spi- 
ders while increasing their number is greater than in the 
case of a classical coronagraph. When designing ELT 
secondary support structures and planning to correct for 
them using ACAD, increasing the number of spiders to 
8 or even 12 has a beneficial impact on contrast as it en- 
ables each discontinuity to become thinner and thus to be 
corrected to higher contrast using the DMs. The PSFs of 
apertures with more than 12 spiders present diffraction 
structures which are poorly suited for correction with 
square DMs. While the contrast resulting from applying 
ACAD to such apertures is still a decreasing function of 
the number of spiders, Fig. [18] shows that the net con- 
trast gain brought by the DIVLbased remapping is smaller 
than in the more gen tle cas es o f 12 spiders. The study 
presented on Fig. [15] to Fig. [18] remains to be fully opti- 
mized for each potential design of an ELT planet finding 
instrument (in particular using a finer grid of cutoff spa- 
tial frequencies). It however demonstrates the flexibility 
of ACAD for various aperture geometries and provides 
a first order rule of thumb to design telescope apertures 
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Figure 15. PSFs resulting from ACAD when varying the number and thickness of secondary support structures while maint aini ng their 
covered surface constant. The surface area covered in this example is 50% greater than in the TMT example shown on Fig. |13| As the 
spiders get thinner their impact on raw contrast becomes smaller and the starlight suppression after DM correction becomes bigger. For 
a relatively small number of spiders (< 12) the contrast improvement on each single structure is the dominant phenomenon, regardless of 
the number of spiders. ELTs designed with a moderate to large number of thin secondary support structure (6 to 12) present aperture 
discontinuities which are easy to correct with ACAD. 
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Figure 16. Radial PSF profiles resulting from ACAD when vary- 
ing the number and thickness of secondary support structures while 
maintaining their covered surface constant. The surface area cov- 
ered in this e xam ple is 50% greater than in the TMT example 
shown on Fig. |13| As the spiders get thinner their impact on raw 
contrast becomes lesser and the starlight suppression after DM 
correction becomes greater. In the 12 spiders example, at large 
separations, the average contrast is an order of magnitude higher 
than reported on Fig. 1 1 4| 



which are friendly to direct imaging of exo-planets: "A 
lot of thin spiders is more favorable than a few thick spi- 
ders" . In practice the number of spiders will be limited 
by effects not treated in our analysis such as the me- 



chanical rigidity, requirements on the perfection of their 
periodic spacing and glancing reflections from the sides 
of multiple spiders. We thus advocate that, should future 
ELTs be built with high contrast exo-planetary science 
as a main scientific driver, then such effect ought to be 
thoroughly analyzed as a large numbers of thin spiders 
is more favorable from a contrast standpoint when using 
ACAD to mitigate for pupil amplitude asymmetries. 

6.2.2. Monolithic on-axis apertures. 

When discussing the case of JWST we stressed the 
complexity associated with the optimization of ACAD in 
the presence of aperture discontinuities of varying width. 
Carrying out such an exercise would be extremely valu- 
able to study the feasibility of the direct imaging of exo- 
earth with an on-axis se gmented future flagsh ip observa- 



tory such as ATLAST (postman et al. 2010). However 



such an effort is computationally heavy and thus beyond 
the scope of the present paper, which focuses on intro- 
ducing the ACAD methodology and illustrating using 
key basic examples. 

So far, none of the examples in this manuscript demon- 
strate that ACAD can yield corrections all the way down 
to the theoretical contrast floor that is set by the coron- 
agraph design. When seeking to image exo-earths from 
space, future missions will need to reach this limit. In or- 
der to explore this regime, we conducted a detailed study 
of an hypothetical on-axis monolithic telescope with four 
secondary support struts. To establish the true contrast 
limits we varied the thickness of the spiders and for each 
geometries. The surfaces were calculated by solving the 
reverse problem over an even grid of 70 cutoff low-pass 
spatial frequencies ranging between 30 and 70 cycles per 
apertures for the tapering kernel. The value yielding the 
best contrast was chosen. In all cases we used a corona- 
graph with a central obscuration of 0.15 D, with a focal 
plane mask of 6 X/D diameter, and IWA of 4X/D, and 
an OWA of 30 X/D. Note that when using coronagraphs 
relying on pupil apodization these results can be readily 
generalized to larger circular secondary obscurations, at 
a loss in IWA (as shown on Fig. [2j). Moreover we clocked 
the telescope aperture by 45° with respect to the grid 
of Fourier modes. We found that, indeed, the theoreti- 
cal contrast floor set by the coronagraph design is met 
for thin spiders (0.02 D) and it is very close to be met 
for spiders only half the thickness of the the ones cur- 
rently equipping the Hubble Space Telecope (0.05 D), 
see Fig. [19] and Fig. [20j Even in the case of thick struts 
(0.1 D) we find contrasts an order of magnitude higher 
tha n in the similar configuration on the top panel of 
Fig. [l6j due to our thorough optimization of the cutoff 
frequency of the tapering kernel and careful clocking of 
the aperture with respect to the actuators. On-axis tele- 
scopes are thus a viable option to image earth-analogs 
from space: their secondary support structures can be 
corrected down to contrast levels comparable to the tar- 
get contrast of recent mission s concept studies (Guyon 
et al.|2008|[Trauger et al.|2010l) . Since the baseline wave- 
front control architecture tor future space coronagraphs 
relies on two sequential DMs, ACAD does not add any 
extra complexity to such missions and merely consists of 
controlling the DMs in order to optimally compensate 
for the effects of asymmetric aperture discontinuities. 
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Figure 17. PSFs resulting from ACAD when varying the number and thickness of secondary support structures while ma inta ining their 
covered surface constant. The surface area covered in this example is 50% greater than in the TMT example shown on Fig. |13| When the 
number of spiders increases, they produce a sharp circular diffraction feature at Ng p i^ ers /TT X/D. If this number is greater than the size 
of the focal plane mask this structure appears in the high contrast zone and is very difficult to correct with ACAD. The brightness of this 
structure is mitigated by the fact that the spiders are very thin. 



24 



Pueyo & Norman 



A/D in case of a PIAAC 
8.7 13. 



DMs at rest 




Correction with N = 64 


' % \ V / 




Ns = 16 






width =0.3 % 




\ A / 





PIAA coronagraph ( Martinache et al. 2006). However 
the asphericity of the surfaces in the case oi ACAD oper- 
ating on reasonably thin discontinuitie s, is much smaller 
than in a PIAA remapping unit. Fig. |21| shows the im- 
pact on off-axis PSFs of such a distortion in the worse- 
case scenario of a geometry similar to JWST. We demon- 
strate that most of the flux remains in the central disk of 
radius X/D for all sources in the field of view of the coro- 
nagraphs considered here (all the way to 30 cycles per 
capture). We conclude that, because of the small defor- 
mations of the DMs, PSF distortion will not be a major 
hindrance in exo-planet imaging instruments whose DMs 
are controlled in order to mitigate for discontinuities in 
the aperture. 
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Figure 18. Radial PSF profiles resulting from ACAD when vary- 
ing the number and thickness of secondary support structures while 
maintaining their covered surface constant. The surface area cov- 
ered in this e xam ple is 50% greater than in the TMT example 
shown on Fig. |13| With a large number of spiders the bright ring 
in the PSF structure located at AfspWers/ 71 " A/D is difficult to cor- 
rect with ACAD. However since the spiders becomes thinner its 
net effect on contrast after ACAD remains small. 



7. DISCUSSION AND FUTURE WORK 

7.1. Field dependent distortion 

Because ACAD relies on deforming the DMs surfaces 
in an aspherical fashion, off-axis wavefronts seen through 
the two DMs apparatus will be distorted, just as in a 



7.2. Impact of wavefront discontinuities in segmented 
telescopes. 

7.2.1. General equations in the presence of incident phase 
errors and discontinuities 



So far we have treated primary mirrors' segmentation 
as a pure amplitude effect. In reality the contrast 
floor in segmented telescopes will be driven by both 
phase and amplitude discontinuities: here we explore 
the impact of phase errors and discontinuities occur- 
ring before two DMs whose surfaces have been set 
using ACAD. There are two main phenomena to be 
considered. The first is the conversion of the incident 
wavefront phase before DM1: ^Ahi into amplitude at 
the second mirror. The second is the projection of this 
wavefront phase into a remapped phase errors at DM2: 
^-Ah 1 (f 1 (x 1 ,y 2 ),gi{x 1 ,y 2 ))- Since the remapping unit 
is designed using deformable mirrors, both DM1 and 
DM2 a complete correction could be attained in prin- 
ciple. However, the deformable mirrors are continuous 
while Ahi presented discontinuities. Thus, complete 
corrections for segmented mirrors might not be achieved 
in practice. Below we discuss the following two main 
points. (1) Even if the phase wavefront error Ahi has 
discontinuities, the phase errors within in segment still 
drive the phase to amplitude conversion and thus the 
propagated amplitude at DM2. In that case treatments 
of these phenomenons that have already been discussed 
in the literature for monolithic aper tures a re still valid 
for small enough phase errors Eqs. |56[56| and smooth 



enough remapping functions. For ACAD remapping 
this smoothness constraint is naturally enforced by the 
limited number of actuators across the DM surface. In 
this case phase to amplitude conversion between can in 
principle be corrected using DM1. (2) Remapped phase 
discontinuities can be corrected for a finite number 
of spatial frequencies using a continuous phase sheet 
deformable mirror. We illustrate this partial correction 
over a 20% bandwidth using numerical simulations of 
a post-ACAD half dark hole created by superposing a 
small perturbation, computed using a linear wavefront 
control algorithm, to the ACAD DM2 surface. 
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Figure 19. PSFs resulting from ACAD when varying the number and thickness of secondary support structures. As the spiders get 
thinner their impact on raw contrast becomes lesser and the starlight suppression after DM correction becomes greater. In this case 10 was 
optimized on a very fine grid and the aperture we clocked in a favorable direction with respect to the Fourier basis. 
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If the incoming wavefront is written as A/ii and the 
solution of the Monge Ampere Equation for DM1 as h\ 
then one can conduct the analysis in Eq. [6] to 12b using 
hi = h® + Ah\. Under the assumption that surface of 
DM2 is set as hi in order to conserve Optical Path Length 
then one can re- write the remapping as (/i,<fi) defined 
by: 



dhi 
dx 
dhi 
dy 



{h(x2,V2),gi(x2,V2)) Z 

{h{^2,V2),gi(x2m)) Z 



(56) 
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Figure 20. PFSs resulting from ACAD when varying the num- 
ber and thickness of secondary support structures. As the spi- 
ders get thinner their impact on raw contrast becomes lesser and 
the starlight suppression after DM correction becomes greater. In 
this case ui was optimized on a very fine grid and the aperture we 
clocked in a favorable direction with respect to the Fourier basis. 
Even for spiders as thick as 0.5% of the telescope aperture the 
designed contrast of the coronagraph is retrieved. 





Figure 21. Off-axis PSF after ACAD in the case of a geometry 
similar to JWST. The aspheric surface of the DMs introduce a 
slight field-dependent distortion. However the core of the PSF is 
still concentrated within the central airy disk and the DMs only 
have an effect on the PSF tail. Field distortion does not thus 
hamper the detectability of faint off-axis sources. 



J 



Moreover if edge ringing has been properly mitigated then the ray optics solution is valid and the field at DM2 can 
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be written as: 



E D M2(X2,V2) = 




E 



DM1 



d 2 hi \ _ / d 2 hx \2 
dy 2 * ^ dxdx * 



J^(S(fi ,ffl )+hi (/i ,g\)-h2 ) 



(/l,9l) 



(57) 



7.2.2. Impact on the amplitude after ACAD 

We first consider the amplitude profile in Eq. [57j it 
is composed of two factors the remapped telescope aper- 
ture, Er>Mi(fi, <7i), and the determinant of /d + 'H^i]. 

The first condition necessary for the incoming wave- 
front not to perturb the ACAD solution is: Ahi is such 
that the remapping is not modified at the pupil locations 
where the telescope aperture is not zero Edmi ^ 0. This 
results into the conditions 

dAhi dAhl , s , „ 

<C — — for (a;,y)sucn thatEnMi {x, y) f= 



dx dx 
dAhi dAh\ 



dy 



dy 



for (x, y)such that^cMi (x, y) ^ 



At the locations where Edmi — there is no light 
illuminating the discontinuous wavefront and thus the 
large local slopes at these location have no impact on 
the remapping functions (/i,<?i). These conditions are 



not true in segmented telescopes that are not properly 
phased, for which the tip-tilt error over each segment 
can reach several waves. However under the assump- 
tion that the primary has been properly phased (for in- 
stance the residual rms wavefront after phasing is expect 
to be ~ 1/10 th of a wave, similar to values expected for 
JWST NIRCAM) then these conditions are true within 
the boundaries of each segment. Moreover, while the lo- 
cal wavefront slopes at the segment's discontinuities do 
not respect this condition the incident amplitude at these 
points is EoMi{x,y) — and they thus do not perturb 
the ACAD remapping solution. 

The second necessary condition resides in the fact that 
the determinant of Id + T-L(h\) is not equal to det[Id + 
H(hi)] at the pupil locations where the telescope aperture 
is not zero ought not have a s evere impact on contrast . 
One can use the linearization in Loeper & Rapetti ( 2005 1 
to show that: 



det[Id + H(h\)} det[Id + H(h1)](l 
1 1 



1 

(l+''L)^''i m +(l+h; !)s )AJ il; , I -2«h; il) A/ illa 
det[Id+H(h1)] 



det[Id + H(hi)] det[Id + H{h\)]{l + AA(Ahi)) 



(58) 



(59) 



The perturbation term AA{Ah±) corresponds to the 
full non-linear expression of the phase to amplitude con- 
version of wav efront errors t hat o ccurs in pupil remap- 
ping units. In |Pueyo et al. (2011) we derived a similar 
expression in the linear case, when Ah\ <C A and showed 
that in the pupil regions where the the beam is con- 
verging this phase to amplitude conversion was enhanced 
when compared t o the case of a Fr esnel propagation. 
In a recent study Krist et al. (2011) presented simula- 
tions predicting that this effect was quite severe in PIAA 
coronagraphs and can limit the broadband contrast af- 
ter wavefront control unless DMs where placed before 
the remapping unit. In principle ACAD will not suffer 
from this limitation as the first aspherical surface of the 
remapping unit is actually a Deformable Mirror that can 
actually compensate for Ahi, before any phase to am- 
plitude wavefront modulation occurs. Devising a wave- 
front controller that relies on DM1 requires moreover a 
computationally efficient model to propagate arbitrary 
wavefronts though t ACAD. Such a tool was developed in 
Krist et al. (20101 assuming azimuthally symmetric ge- 
ometries. Since devising such a tool in ACAD's case, in 
the asymmetric case, represents a substantial effort well 
beyond the problem of prescribing ACAD DM shapes, 
we chose not to include such simulations in the present 
manuscript. Since we are using Deformable Mirrors with 



a limited number of actuators, ACAD remapping is in 
general less severe than in the case of PIAA. We thus 
expect the results regarding the wave front correct i on be- 
fore the remapping unit reported in |Krist et al. (20111 
to hold. This is provided that the DM actuators can 
adequately capture the high spatial frequency content of 
Ahi to create a dark hole in the coronagraphic PSF. We 
tackle this particular aspect next when discussing the 
case of phase errors, in the absence of wavefront phase 
to amplitude conversion. Once again note that while 
the local wavefront curvatures are very large at the seg- 
ment's discontinuities, the incident amplitude at these 
points is EjjMl(x,y) = and they thus they do not 
have an impact on the ACAD phase to amplitude mod- 
ulation. In practice if the DM is not exactly located at 
a location conjugate to the telescope pupil the actual 
wavefront discontinuities will be slightly illuminated and 
might perturb the remapping functions and the phase to 
amplitude conversion. While this might tighten require- 
ments regarding the positioning of DM1 in the direction 
of the optical axis we do not expect this effect to be a 
major obstacle to successful ACAD implementations. 

7.2.3. Impact on the phase after ACAD 

In practice, when Ahi presents discontinuities, the sur- 
face of DM2 cannot be set to the deformation /12 that 
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Figure 22. Broadband wavefront correction (20% bandwidth around 700 nm) with a single DM in segmented telescope with discontinuous 
surface errors. Top Left: wavefront before the coronagraph. Top Right: broadband aberrated PSF with DM at rest. Bottom Left: DM 
surface resulting from the wavefront control algorithm. Bottom Right: broadband corrected PSF. Note that the wavefront control algorithm 
seeks to compensate for the diffractive artifacts associated with the secondary support structures: it attenuates them on the right side of 
the PSF while it strengthens them on the left side of the PSF. As a result the DM surface becomes too large at the pupil spider's location 
and the quasi-linear wavefront control algorithm eventually diverges. 
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conserves Optical Path Length, since we work under the 
assumption that the DMs has a continuous phase-sheet. 
While this has no impact on the discussion above regard- 
ing the amplitude of Edm2, since DM1 is solely responsi- 
ble for this part, it ought to be taken into account when 
discussing the phase at DM2. Under the assumption 
that Ah\ does not perturb the nominal ACAD remap- 
ping function then one can show that the phase at DM2 
is: 

arg[E DM2 (x 2 ,y 2 )} = 
= —(Ah 1 (f^(x 2 ,y 2 ), s°(* a , y 2 )) + Ah 2 {x 2 , y 2 )) (60) 

where Ah 2 (x 2 ,y 2 ) is a small continuous surface de- 
formation superposed to the ACAD shape of DM2 and 
Ah 1 (f^(x 2 ,y 2 ),g^(x 2 ,y 2 )) is the telescope OPD seen 
through the DM based remapping unit. This second 
term presents phase discontinuities whose spatial scale 
has been contracted by ACAD. When these discontinu- 
ities are very small then their high spatial frequency con- 
tent does not disrupt the ability of DM2 to correct for low 
to mid-spatial frequency wavefront errors wavefront er- 
rors. However as the discontinuities become larger their 
high spatial frequency content can fold into the region 
of the PSF that the DMs seek to cancel. These "fre- 



quency fold ing" speckles are highly chromatic (Give'on 
et al. 2006) and can have a severe impact on the spectral 



bandwidth of a coronagraph whose wavefront is corrected 
using a continuous DM. 

In order to assess the impact of this phenomenon, we 
conducted a scries of simulations based on single DM 
wavefront control algorithm that seeks to crea te a dark 
hole in one h alf of the image plane at in as in |Borde fc] 
Traub ( 2006[) . We use the example of a geometry sim- 
liar to J WST and work under the assumption that the 
discontinuous wavefront incident to the coronagraph has 
the same spatial frequency content as a JWST NIRCAM 
Optical Path Difference that has been adjusted to 70 
nm rms in order to mimic a visible Strehl similar to the 
near-infrared Strehl of JWST. The non-linear wavefront 
and sensing and control problem associated with phasing 
a primary mirror to such level of precision is undoubt- 
edly a colossal endeavor and is well beyond the scope of 
this paper. In this section we work under the assump- 
tion that the primary mirror either has been phased to 
such a level, that the wavefront discontinuities are no 
larger than 200 nm peak to valley or that the wavefront 
has been otherwise corrected down to this specification 
using a segmented Deformable Mirror that is conjugate 
with the primary mirror. Moreover we assume (1) that 
the residual post-phasing wavefront map has been char- 
acterized and can be used in order to build the linear 
model underlying the wavefront controller (2) the focal 
plane wavefront estimator ( carried using DM diversity 
as in Borde & Traub (2006) for instance) is capable to 
yield an exact estimate of the complex electrical field at 
the science camera. Underlying this last assumption is 
the overly optimistic premise that wavefront will remain 
unchanged over the course of each high-contrast expo- 
sure. While this is not a realistic assumption one could 
envision the introduction of specific wavefront sensing 
schemes, with architectures similar to the one currently 
considered for low order wavefront sensors on monolithic 



apertures ( |Guyon et al.|2009| |Wallace et al.|2011[ ) , or us- 
ing a separate metrology system. The results presented 
here are thus limited to configurations for which segment 
phasing will be dynamically compensated using specific 
sensing and control beyond the scope of this paper. As 
this section merely seeks to address the controllability of 
wavefront errors in segmented telescopes we chose to con- 
duct our simulations with a perfect estimator. Finally we 
use the strok e minimization wave front control algorithm 
presented in Pueyo et al. (2009) to ensure convergence 
for as many iteration as possible. We first tested the case 
of a segmented telescope in the absence of ACAD, using 
a azimuthally symmetric coronagraph and a single DM. 
We sough to create a Dark Hole between 5 and 28 Xq/D 
under a 20% bandwidth with A 



700 nm. Fig 22 



shows 

the results of such a simulation. The DM can indeed 
correct for the discontinuities over a broadband in one 
half of the image plane. However the wavefront control 
algorithm seeks to compensate for the diffractive arti- 
facts associated with the secondary support structures: 
it attenuates them on the right side of the PSF while it 
strengthens them on the left side of the PSF. As a result 
the DM surface becomes too large at the pupil spider's lo- 
cation and the quasi-linear wavefront control algorithm 
eventually diverges for contrasts ~ 10 6 . We then pro- 
ceeded to simulate the same configuration in the pres- 
ence of two DMs whose surface at rest was calculated 
using ACAD. Since there does not exist a model yet 
to propag ate arbitrary wave fronts through ACAD (the 



models in Krist et al. ( 2011[ ) only operate under the as- 
sumption of an azimuthally symmetric remapping) we 
can only use the second DM for wavefront control. We 
work under the assumption that the incident wavefront 
does not perturb the nominal ACAD remapping (which 
is true in the case of the surface ma p we chose f or our 



example) and that the arguments in |Krist et al. ([20101 
hold so that phase to amplitude conversion m ACAD 
can be compensated by actuating DM1. Frequency fold- 
ing will then be the phenomenon responsible for the true 
contrast limit. In this section we are interested in ex- 
ploring how this impacts the controllability of wavefront 
discontinuities using continuous phase-sheet DMs. We 
used a azimuthally symmetric coronagraph and super- 
posed our wavefront control solution to DM2. We sough 
to create a Dark Hole between 5 and 28 Xo/D under a 
20% bandwidth with Ao = 700 nm. Fig 22 shows the re- 
sults of such a simulation. When the incident wavefront 
is small enough it is indeed possible to superpose a "clas- 
sical linear wavefront control" solution to the non-linear 
ACAD DM shapes in order to carve PSF dark holes. 
The wavefront control algorithm now yields a DM sur- 
face that does not feature prominent deformations at the 
location of the spiders. Most of the DM stroke is located 
at the edge of the segments, at location of the wave- 
front discontinuities and seek to correct the frequency 
folding terms associated with such discontinuities. At 
these locations the DM surface eventually becomes too 
large and the linear wavefront control algorithm diverges. 
However this divergence occurs at contrast levels much 
higher than when the ACAD solution is not applied to 
the DMs. These simulations show that indeed discon- 
tinuous phases can be corrected using the second DM 
of a ACAD whose surfaces have preliminary been set to 
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Figure 23. Broadband wavefront correction (20% bandwidth around 700 nm) in a segmented telescope whose pupil has been re-arranged 
using ACAD. The surface of the first DM is set according to the ACAD equations. The surface of the second DM is the sum of the 
ACAD solution and a small perturbation calculated using a quasi-linear wavefront control algorithm . Top Left: wavefront before the 
coronagraph. Note that the ACAD remapping has compressed the wavefront errors near the struts and the segment gaps. Top Right: 
broadband aberrated PSF with DMs set to the ACAD solution. Bottom Left: perturbation of DM2's surface resulting from the wavefront 
control algorithm. Bottom Right: broadband corrected PSF. The wavefront control algorithm now yields a DM surface that does not 
feature prominent deformations at the location of the spiders. Most of the DM stroke is located at the edge of the segments, at location 
of the wavefront discontinuities. There, the DM surface eventually becomes too large and the quasi-linear wavefront control algorithm 
diverges. However this occurs higher contrasts than in the absence of ACAD. 



mitigate the effects of spiders and segment gaps. 

7.3. Ultimate contrast limits 

Assuming that edge ringing has been properly miti- 
gated, so that the ray optics approximation underlying 
the calculation of the DMs shapes is valid, one can won- 
der about the ultimate contrast limitations of the results 
presented in this manuscript. Increasing the number of 
actuators would have dramatic effects on contrast if the 
actuator count would be such that N > D/d where d is 
the scale of the aperture discontinuities. Unfortunately 
current DM technologies are currently far from such a 
requirement and the solutions presented here are in the 
regime where N <C D/d. In this regime N only has 
a marginal influence on contrast when compared to the 
impact of the cutoff frequency of the tapering kernel. 



In the regime described here varying the actuator count 
only changes the size of the corrected region. 

The residual PSF artifacts in Figs, [ft] to [20] follow 
the direction of the initial diffraction pattern associated 
with secondary support structures and segments. When 
addressing the problem of aperture discontinuities by 
solving the Monge-Ampere Equation, ACAD calculates 
the DM shapes based on a pupil plane metric and thus 
mostly focuses on attenuating these structures with 
little regard to the final contrast. It is actually quite 
remarkable that such a pupil-only approach yields levels 
of starlight extinction of two to three orders of mag- 
nitude. A more appropriate metric would be the final 
intensity distribution in the post-coronagraphic image 
plane. However, as discussed in §. [3] classical wavefront 
control algorithms based on a linearization of the DMs 



High contrast with arbitrary apertures 



31 



5.8 



X/D in case of a PIAAC 
11.5 17.2 



Aberrated ACAD PSF 
Corrected ACAD PSF 
ACAD PSF with no aberratio 



Aberrated non ACAD PSF 
Corrected non ACAD PSF 




Feedback/DM 
control strategy 



System 



Forward 
modeling tools 



X/D in case of an APLC 



Fig ure 24. Ra dial average in the half dark plane of the PSFs on 
Fig |22| and Fig. |23| In the presence of wavefront discontinuities 
corrected using a continuous membrane DM, ACAD still yields, 
over a 20% bandwidth around 700 nm, PSF with a contrast 100 
larger than in a classical segmented telescope. Moreover this fig- 
ure illustrates that since it is based on a true image plane metric, 
the wavefront control algorithm can be used ( within the limits of 
its linear regime) to improve upon the ACAD DM shapes derived 
solving the Monge Ampere Equation. 



deformations around l ocal equilibrium s hapes (such as 
the ones pr esented in Borde fc Traub fl2006|); Give'on 



et al. (2007) in the one DM case or Pueyo et ai 
tor one or two DMs) cannot be used to compensate 
the full aperture discontinuities. This is illustrated in 
Fig [22] where the DM surface in the vicinity of spiders 
becomes too large after a certain number of iterations, 
which leads the iterative algorithm to diverge. When 
attempting to circumvent this problem by recomputing 
the linearization at each iteration, we managed to 
somewhat stabilize the problem for a few iterations and 
reached marginal contrast improvements, but the overall 
algorithm remained unstable unless a prohibitively 
small step size was used. This is the problem which 
motivated our effort to calculate the DM shapes as the 
full non-linear solution of the Monge- Ampere Equation. 
While doing so yields significant contrast improvements 
in both the case of JWST like geometries, TMT and on 
axis-monolithic apertures similar, this approach does 
not give a proper weight to the spatial frequencies of 
interest for high contrast imaging. We mitigated this 
effect by giving a strong weight to the spatial frequencies 
of interest (in the Dark Hole) when solving the Monge 
Ampere Equation. 

The next natural step is thus to use non-linear solu- 
tions presented herein to correct for the bulk of the aper- 
ture discontinuities and to serve as a starting point for 
classical linearized waveform control algorithms, as illus- 
trated on Fig. [25] Fig. [24] indeed illustrates that when 
superposing an image plane based wavefront controller 
to the Monge Ampere ACAD solution, the contrast can 
be improved beyond the floor shown on Fig. |12| How- 
ever one DM solutions, are of limited interest as they 
only operate efficiently over a finite bandwidth and over 
half of the image plane. ACAD yields a true broad- 
band solution, and consequently it would be preferable 
to use the two DMs in the quasi-linear regime to quan- 




Figure 25. Future work towards higher contrasts with ACAD. 
The blue and orange colors respectively represent the current state 
of the art in wavefront control and the work described in the present 
manuscript, as in Fig.JT] In brown are listed the potential avenues 
to further the contrasts presented herein: 1) combining ACAD with 
coronagraphs designed on segmented and/or on-axis apertures, 2) 
using diffractive models to close a quasi-linear focal plane based 
loop using a metric whose starting point corresponds to the DM 
shapes calculated in the non-linear regime. 

tify the true contrast limits of ACAD . In such a scheme 
the DM surfaces are first evaluated as the solution of 
the Monge-Ampere Equation and then adjusted using 
the image plane based wavefr ont control algorithm pre- 
sented in |Pueyo et al. ( 2009 1 . However such an exer- 
cise requires effic ient and robust numerical algorithms to 
evaluate Eq. [16} Such tool only exist so far in the case 
of az imuthally symmetric remapping units (Krist et al. 
2010 1. Developing such numerical tools is thus ot pri- 



mary interest to both quantifying the chromaticity and 
the true contrast limits achievable with on-axis and/or 
segmented telescopes. 

8. CONCLUSION 

We have introduced a technique that takes advantage 
of the presence of Deformable Mirrors in modern high- 
contrast coronagraph to compensate for amplitude dis- 
continuities in on-axis and/or segmented telescopes. Our 
calculations predict that this high throughput class of so- 
lutions operates under broadband illumination even in the 
presence of reasonably small wavefront errors and discon- 
tinuities. Our approach relies on controlling two sequen- 
tial Deformable Mirrors in a non-linear regime yet unex- 
plored in the field of high-contrast imaging. Indeed the 
mirror's shapes are calculated as the solution of the two- 
dimensional pupil remapping problem, which can be ex- 
pressed as a non-linear partial differential equation called 
Monge Ampere Equation. We called this technique Ac- 
tive Compensation of Aperture Discontinuities. While 
we illustrated the efficiency of ACAD using Apodized 
Pupil Lyot and Phase Induced Amplitude Coronagraph, 
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it is is applicable to all types of coronagraphs and thus 
enables one to translate the past decade of investigation 
in coronagraphy with unobscured monolithic apertures 
to a much wider class of telescope architectures. Be- 
cause ACAD consists of a simple remapping of the tele- 
scope pupil, it is a true broadband solution. Provided 
that the coronagraph chosen operates under a broadband 
illumination, ACAD allows high contrast observations 
over a large spectral bandwidth as pupil remapping is 
an achromatic phenomenon. We showed that wavelength 
edge diffraction artifacts, which are the source of spec- 
tral bandwidth limits in PIAA coronagraphs (also based 
on pupil remapping), are no larger than classical Fresnel 
ringing. We thus argued that they will only marginally 
impact the spectral bandwidth of a coronagraph whose 
input beam has been corrected with ACAD. The mirror 
deformations we find can be achieved, both in curvature 
and in stroke, with technologies currently used in Ex- AO 
ground based instruments and in various testbeds aimed 
at demonstrating high-contrast for space based applica- 
tions. Implementing ACAD on a given on-axis and/or 
segmented thus does not require substantial technology 
development of critical components. 

For geometries analogous to JWST we have demon- 
strated that ACAD can achieve at least contrast ~ 10~ 7 , 
provided that dynamic high precision segment phasing 
can be achieved. For TMT and ELT, ACAD can achieve 
at least contrasts ~ 1CP 8 . For on-axis monolithic obser- 
vatories the design contrast of the coronagraph can be 
reached with ACAD when the secondary support struc- 
tures are 5 times thinner than on HST. When they are 
just as thick as HST contrasts as high as 10 8 can be 
reached. These numbers are, however, conservative: an 
optimal solution can be obtained by fine tuning the con- 
trol term in the Monge Ampere Equation to the charac- 
teristic scale of each discontinuity. As our goal was to 
introduce this technique to the astronomical community 
and emphasize its broad appeal to a wide class of ar- 
chitectures (JWST,ATLAST,HST,TMT,E-ELT) we left 
this observatory specific exercise for future work. 

The true contrast limitation of ACAD resides in the 
fact that the Deformable Mirrors are controlled using a 
pre-coronagraph pupil based metric. However, as illus- 
trated in Fig. 25 the solution provided by ACAD can be 
used as the starting point for classical linearized wave- 
form control algorithms based in image plane diagnostics. 
In such a control strategy, the surfaces are first evalu- 
ated as the solution of the Monge- Ampere Equation and 
then adjusted using the quasi-linear method presented in 
Pueyo et al.| ([2009 ) . This control strategy requires effi- 
cient and robust numerical algorithms to evaluate the full 
diffractive propagation in the remapped Fresnel regime. 
All the contrasts reported here are achieved without 
aberrations and we showed that in practice, quasi-linear 
DM controls based on images at the science camera will 
have to be superposed to the ACAD solutions. Finally, 
as ACAD is broadly applicable to all types of corona- 
graphs, the remapped pupil can be used as the entry 
point to relax the design of coronagraphs that do oper- 
ate on segme n ted apertures such as discussed in |CarlottI| 



et al. (2011); 



Fig. p5 



Guyon et al. (2010a), also illustrated in 



put broadband solutions for a variety of coronagraphs. It 
only relies on hardware (Deformable Mirrors) that have 
been extensively tested over the past ten years. Finally 
since ACAD can operate with all type of coronagraphs 
and it renders the last decade of research on high-contrast 
imaging technologies with off-axis unobscured apertures 
applicable to much broader range of telescope architec- 
tures. 
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ACAD is thus a promising tool for future high 



contrast imaging instruments on a wide range of observa- 
tories as it will allow astronomers to devise high through- 
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